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Abstract 



Charge conjugation C is one of the fundamental symmetries in nature which transforms 
particles into antiparticles. This symmetry was studied in weak interaction where it is fully 
violated, but it is poorly known in the strong and electromagnetic interactions. Therefore, 
it is important to test this symmetry accurately for a better understanding of the nature 
of the strong interaction and for the understanding of the significantly larger abundance 
of matter over antimatter in the Universe. 

To this end, in this thesis we investigated r] — >■ Tr+Tr^Tr'^ and 77 — >■ Tr^e+e" decays, which 
might violate charge conjugation symmetry. The violation of C symmetry in 77 — >■ 7r"'"7r~7r^ 
process could manifest itself as an asymmetry between energy spectra of charged pions, 
and can be studied using event density distribution on Dalitz plot. The 77 — >■ Tr^e'^e~ 
decay is forbidden by C symmetry in the first order of electromagnetic interaction, and 
can only proceed by emission of two virtual photons with the branching ratio on a level 
of 10~^, as predicted in the framework of the Standard Model. Therefore, observation of 
a larger branching ratio could indicate a mechanism involving first order electromagnetic 
interaction violating charge conjugation. 

Both decays were investigated by means of the WASA-at-COSY detector operating 
at the COSY synchrotron at the Forchungszentrum Jiilich in Germany. The r] meson 
was produced via pp — pprj reaction at the proton beam momentum of 2.14 GeV/c which 
corresponds to kinetic energy of 1.4 GeV. The measurement was done at the turn of October 
and November in the year 2008. In total around 5 • 10^ rj mesons were collected during the 
two weeks of data taking. The tagging of the 77 meson was done by means of the missing 
mass technique and the decay products were identified by the invariant mass reconstruction. 

As a result of the analysis conducted in the framework of this thesis a Dalitz Plot 
distribution for the t] — )■ tt~^'k~7:^ decay was obtained. From this distribution we extracted 
asymmetry parameters sensitive to C symmetry violation for different isospin values of the 
final state and we have established that all are consistent with zero within the obtained 
accuracy. 

For the r) 7r°e"'"e~ decay we have not observe signal and thus we estimated an upper 
limit for the branching ratio. The established upper limit amounts to BRj^^^q^+^- < 
3.7 • 10~^ at the 90% confidence level. This result is more precise than previously obtained 
in other experiments. We intend to continue the research, and thanks to the 20 times higher 
statistics of already collected data by WASA-at-COSY, the upper limit will be improved 
significantly. 



Streszczenie 



Sprz^zenie ladunkowe C jest jedn^ z podstawowych symetrii w przyrodzie, ktora za- 
mienia czqstki na antycz^stki. Symetria ta byla badana w oddzialywaniach slabych w 
ktorych jest calkowicie lamana, natomiast do tej pory jest slabo poznana z punktu widzenia 
oddzialywan silnych i elektromagnetycznych. Dlatego waznym aspektem jest badanie stop- 
nia zachowania tej symetrii dla lepszego zrozumienia natury oddzialywania silnego oraz 
wyjasnienia wi^kszej abundancji materii niz antymaterii we Wszechswiecie. 

W tym celu zbadalismy dwa procesy rj — >■ 7r''"7r'~7r° i ry — >■ 7r°e'*"e~, ktore mogq. la- 
mac symetri^ sprz^zenia ladunkowego. Niezachowanie sprz^zenia ladunkowego w proccsic 
r) — >■ 7r"'"7r~7r'^, moze ujawnic si^ jako asymetria pomi^dzy rozkladami energii pionow nalad- 
owanych i moze zostac zaobserwowana za pomoc^ badania g^stosci obsadzcn na wykresie 
Dalitza. Natomiast rozpad 7] — t- vr'^e^e^ jest zabroniony przcz symctriq ladunkow^ w 
pierwszym rz^dzic oddzialywan elektromagnetycznych i moze zachodzic tylko przez emisje 
dwoch wirtualnych fotonow ze stosunkiem rozgal^zieh na poziomie 10~^, wedlug przewidy- 
wan w ramach Modelu Standardowego. Jednkaze, zaobserwowanie wi^kszego stosunku 
rozgal^zien swiadczyloby o innym mechanizmie reakcji niz ten przewidywany na gruncie 
Modelu Standardowego, ktory moglby niezachowywac symetrii sprz^zenia ladunkowego. 

Oba rozpady byly badane detektorem WASA-at-COSY zainstalowanym na sychrotronie 
COSY znajdujqcym si^ w Centrum Badawczym Jiilicli w Niemczech. Mezon 77 byl pro- 
dukowany w reakcji pp — t- ppr] przy p^dzie wi^zki protonowej 2.14 GeV/c, co odpowiadalo 
energii kinetycznej 1.4 GeV. Pomiar zostal wykonany na przelomie pazdziernika i listopada 
2008 roku. W trakcie dwoch tygodni zebrano probk^ danych zawieraj<^c<^ okolo 5 • lO'^ 
mezonow rj. Mezon r] byl identifikowany za pomoc^ widma masy brakuj^cej, natomiast 
produkty jego rozpadu byly zidentyfikowane za pomoc^ rekonstrukcji masy niezmienniczej . 

W oparciu o przeprowadzonq, w ramach pracy doktorskiej analiz^ danych doswiadczal- 
nych zrekonstruowano wykres Dalitza. Na jego podstawie oszacowano wartosci parametrow 
asymetrii czulych na niezachowanie symetrii sprz^zenia ladunkowego C dla roznych wartosci 
izospinu czq,stek w stanie kohcowym. Stwierdzono, ze wszystkie otrzymane wartosci asymetrii 
sq. zgodne z zerem w granicach oszacowanych niepewnosci. 

Dla rozpadu rj — >■ 7r°e'*"e~ w wyniku przeprowadzonej analizy nie zaobserwowano syg- 
nalu i dlatego oszacowano gornq granic^ stosunku rozgal^zieh na ten rozpad. Obliczona 
wartosc gornej granicy wynosi BR^^^o^+^~ < 3.7 • 10^^ na poziomie ufnosci 90%. Wynik 
ten jest bardziej dokladny niz uzyskany w poprzednich eksperymentach. W najblizszej 
przyszlosci dzi^ki zebranej do tej pory 20 krotnie wi^kszej probce danych przez grup^ 
WASA-at-COSY wynik ten moze zostac znacz^co poprawiony. 
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1. Introduction 



Mesons - the states of the quarl|^and anti-quark (qq) - for more then 60 years play an 
important role in experimental and theoretical physics. Previously and nowadays physicists 
use different mesons to study limits of applicability of the Standard Model [l}|6], which 
is well established theory describing the strong, electromagnetic and weak interactions 
between elementary particles. The examination of mesons production and their decay 
modes give a possibility to probe fundamental symmetries such as: charge conjugation 
(C), space reflection (P), time reversal (T) and their combinations: CP |7| and CPT [Sjjo]. 
Moreover investigation of such processes can be used to determine the parameters of the 
Standard Model. 

One of the particle used for these studies is the rj meson discovered in Berkeley Beva- 
tron Laboratory in 1961 |10| . It is one of the Goldstone bosons in the quantum chromo- 
dynamics (QCD) |11| with no electric charge, flavorless and the mass of rrir] = 547.853 it 
0.024 MeV |12] . From the theoretical point of view it is a superposition of the SU(3) octet 
rjs and singlet rji states, with the wave function: 

= cos 0|r/8) — sin G|7yi), (1-1) 

where G denotes the pseudo-scalar mixing angle between singlet and octet stat^ and: 

{rjs) = —=\uu + dd — 2ss) , (1-2) 

Irji) = —=\uu + dd + ss) . (1-3) 
v3 

The T] meson belongs to the pseudo-scalar family with isospin and angular momentum 
equal to zero, negative parity and charge conjugation equal to +1 {I^{J^'~") = 0''~(0 '')), 
along with the ry', tt^ , vr"*", vr" and , K~, K^, mesons. It is an eigenstate of the 
charge conjugation (C) and parity (P) operators, and thus it constitutes an important 
experimental tool for investigations of the degree of conservation of these symmetries in 
strong and electromagnetic interactions. Moreover, the total width of the r] meson is equal 
to r.^ = 1.30 lb 0.07 keV |12| which is five orders of magnitude smaller than the typical 
width of neutral particles which may decay due to the strong interaction. Therefore, the 
decays of the rj meson are very sensitive to C and P violation. 

The main decay modes of the r] meson can be divided into two groups: hadronic decays 
and radiative decays. All of these strong and electromagnetic processes are forbidden in 



^Quarks are elementary particles with electric charge g = ±| org = ±|. 

^The pseudo-scalar mixing angle was established in [IS and amounts to 6 ~ —15.5° ± 1.3°. 
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the first order |14|. The most energetically favorable strong decays into tt+vt^ and tt^tt^ are 
P and CP violating, with predicted very low branching ratios. Moreover, strong decay into 
47r is also forbidden because of the P and CP invariance and a small available phase space. 
Therefore, the r] meson decays predominantly into tttttt although this decay violates isospin 
symmetry and G-parity 15 . Historically the r] meson decays into tttttt were considered as 



electromagnetic process but it was shown that these effects are small |16| , and instead it is 
expected that the decay occurs only due to the difference in the mass of the u and d quarks. 
This fact permits to study this mass difference by comparing the measured branching ratios 
with the calculations based on the Chiral Perturbation Theory (ChPT) |17|. Furthermore, 
the first order electromagnetic decays like tt^^, tt^tt^^ and 7r+7r^7 are also forbidden and 



they can occur only due to QCD anomalies found in current algebra 18 -201. In the massless 



quark limit the radiative decays are driven by the QCD box anomaly 21 . The second 



order electromagnetic decay r/ — t- 77 is also forbidden, and occurs due to the QCD triangle 
anomaly |22| . The above highlighted properties makes the rj meson specially suitable for 
tests of the discrete symmetries. In this work we use the rj meson for the studies of the C 
symmetry. 

The charge conjugation invariance was studied in weak interaction, and already in 
1957 it was discovered that breaking of this symmetry occurs in decays of vr^ and /j, 



23 



Furthermore, it was also realized that the C operator should turn left-handed neutrinos 
into left-handed anti-neutrinos, but the experimental studies show that all neutrinos are 
left-handed and anti-neutrinos right-handed. This implies that C symmetry should be 
fully violated in weak interactions. It is also interesting to notice that the Big-Bang model 
predicts the same amount of matter and antimatter in the Universe but the experimental 



observations show that there is significantly larger abundance of matter over antimatter 24 



25 . The known CP breaking effect [7| is insufficient to explain this phenomenon, but it 
is hoped that investigations of the charge conjugation may help in clarification of this 
problem. 

One of the main purposes of this thesis is to study the charge conjugation invariance 
in strong interactions by means of the Dalitz plot density population for the r] — )• tt^tt~tt^ 
decay. The r] — )• tt~^tt~tt^ hadronic mode is one of the most frequently occurring decay of 
the rj meson with the branching ratio equal to 22.74 it 0.28% |12j. The C invariance in 
this decay can manifest itself as an asymmetry between the energy distribution of the vr"*" 
and TT~ mesons in the rest frame of the r] meson. The studies of the density population of 
the Dalitz plot can also reveal details of contribution to possible C violation from various 
isospin states of the final particles. Such effects can be investigated by means of three 
asymmetry parameters: (i) A^/j - left-right asymmetry sensitive to violation averaged 
over all isospin states, (ii) Aq - quadrant asymmetry sensitive for the 1 = 2, and (iii) As 



sextant asymmetry which can test the C violation in I = 1 state 26 



Furthermore, we intend to extract the branching ratio or estimate an upper limit for the 
rare rj — ?■ TT^e~^e~ process which might not conserve charge conjugation. In the framework of 
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the Standard Model, the decay rj — )• n^e'^e may only proceed via C-conserving exchange 



of two virtual photons with the branching ratio of about 10~ 27 . But in principle it 
may also be realized by one photon intermediate state, forbidden by the C invariance and 
increasing the branching ratio. At present only an upper limit is set for this branching 



ratio at the level of 4 x 10~^ 12 . Thus, there is still more than three orders of magnitude 
difference between Standard Model predictions and experimentally measured upper limit, 
and therefore increase of the experimental sensitivity gives a chance to observe a signal 
which would indicate violation of C symmetry. The possible charge conjugation breaking 
could be indicated if the branching ratio would be larger than 10^^. 

The measurement aming at the charge conjugation studies described in this dissertation 
was carried out in the Research Center Jiilich in Germany, by means of the WASA-at-COSY 
detector. The r] meson was produced in proton-proton collisions at beam momentum of 
2.14 GeV/c. Identification of the investigated reactions was based on the selection of events 
corresponding to the 7r+7r~7r'^ and Tr^e~^e~ final state. The r] meson signal was extracted 
using missing mass spectrum of two outgoing protons registered in the Forward Detector, 
and the decay products were identified based on invariant mass distribution reconstructed 
from signals detected in the Central Detector of the WASA-at-COSY system. 

In the next chapter the main theoretical motivation for conducted investigations is 
outlined. 

The WASA-at-COSY detector facility and the measurement methods is described in 
Chapter 3. 

The Chapter 4 is devoted to the description of the analysis methods and simulation of 
the detector response. 

In Chapter 5 the track selection methods and reconstruction algorithms are discussed. 

Chapter 6 is committed to the identification method of the pp — )• ppr] reaction. 

Chapter 7 comprises the description of the rj — )■ 7r"'"7r~7r'' decay signal extraction, and 
the discussion of multi-pion background reduction methods. Moreover, the kinematic fit 
procedure is explained in this chapter. 

The final results concerning the Dalitz plot and the asymmetry parameters studies for 
the r] —7- 7r+7r~7r'^ decay are presented in Chapter 8. In addition to that, in Chapter 8 the 
physical background subtraction, the acceptance and efficiency corrections are discussed. 
Finally, the achieved experimental results are compared to theoretical predictions. 

Further on, the procedure of the rj — )• vr'^e^e" decay identification will be presented in 
Chapter 9. 

In Chapter 10, the branching ratio results for the rj — )• vr'^e'^e" decay are discussed. 
The summary and final conclusions followed by perspectives are presented in Chap- 
ter 11. 
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2. Charge conjugation invariance 
tests 



In the Standard Model of particles and fields, the charge conjugation C, along with 
the spatial parity P and the time reversal T, is one of the most fundamental symmetries. 
The C operator in quantum field theory applied to a particle state changes all additive 
quantum numbers of this particle to opposite sign, leaving the mass, momentum and spin 
unchanged, and making it an antiparticle state: 



In the Quantum Electrodynamics (QED) and Quantum Chromodynamics (QCD) it is 
postulated that C holds in all electromagnetic and strong interactions on the level smaller 
than 10~^. Therefore, the C-invariance should imply the balance between matter and 
antimatter, which is not the case in the observed Universe. The Standard Model of the weak 
interactions allows for full C and P violation, as well as for small CP violation. However, 
the model does not explain why and how the violation occurs. Also the discovered small 
CP breaking docs not explain the larger abundance of matter over antimatter. Thus, the 
investigation of the charge conjugation symmetry is one of the most interesting, valuable 
and significant field in modern experimental nuclear and particle physics. 

Difficulties in studies of the charge conjugation arise from the fact that there are only 
very few known particles in nature which are the eigenstates of the C operator. The 
most suitable candidates are neutral and flavorless mesons and the particle-antiparticle 
systems. The particularly interesting appears the 77 meson, which plays a crucial role for 
understanding of the low energy Quantum Chromodynamics, and can be also used to tests 
of the fundamental symmetries. 

In this thesis study of the charge conjugation invariance C is presented, by conducting 
the analysis of the rj — >■ tt'^tt^tt^ and rj — >■ ir^e'^e" decay modes measured in the experiment 
where the rj meson was produced in the proton-proton interactions. 



In view of the tests of charge conjugation invariance C, the strong hadronic decay 
r] — >■ Tr~^7T~7T^ which does not conserve isospin since the Bose symmetry forbids the three 
pions with the = 0~, is particularly interesting. However this decay has to be treated 
in special way because in QCD, at low energies, the strong coupling constant is large. 



civ) = m- 



(2.1) 



2.1 Decay of the r] meson into tt 
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and the perturbative approach is not vahd any more. Therefore, one apphes the Chiral 
Perturbation Theory (ChPT) as an effective field theory specially suited for low energies 
regime. This theory is based on the approximate chiral symmetry and expansion in external 
momenta and quark masses. In this approach the role of dynamic degrees of freedom of 



strong interaction are given to hadrons composed of confined quarks and gluons 28 . One 
can identify them with eight Goldstone bosons members of the pseudoscalar meson nonet, 
which are the result of the spontaneous chiral symmetry breaking [2^. The effective 
Lagrangian is expanded in definite number of derivatives or powers of quark masses given 
as: 

L''JipT = L2 + L4 + Le + ..., (2.2) 

where the subscripts stands for the chiral order. Furthermore, effective Lagrangian shares 
the same symmetries with QCD, namely: C, P, T, Lorentz invariance, and the chiral 
SU{3)l X SU{3)r symmetry. One can see that only even chiral powers arise since the 
Lagrangian is Lorentz scalar which implies that indices of derivatives appear in pairs. The 
lowest order of chiral Lagrangian has only two constants: Bq which is the quark condensate 
parameter, and F^^ which is the pion decay constant, and the decay mechanism is given by 



Current Algebra 30 31] . 

Historically the rj — )• TT~^7r~7r^ decay was treated as an electromagnetic process with 
partial width smaller than second order electromagnetic decay. But as it was shown the 



electromagnetic contributions are small 16 32 and instead the process is dominated by 



the isospin violating term in the strong interaction 33 . Therefore, it is very interesting 



to concern hadronic decays of the j] meson into three pion system in terms of the different 



isospin states. The wave function for 7 = state can be written in the following form 15 



(37r),=o = [(vr+^°),=i|7r-) - (vr+vr-),^!^") + (^-7r°),=i|7r+)] . (2.3) 

where final 3tt can be in the isospin zero state only if the vrvr subsystem is in / = 1 state. 
For the two pion subsystems with 1=1 one can write the wave functions as: 

(vr+^Vi = y|[k+)|vrO)-|vrO)|vr+)], 

(vr+vr-),=i = - IO|vr+)] , 

("°)/=i = /^[-|vr-)k°) + |7rO)|vr-)]. 
Thus the full wave function for the Sir system in / = state reads: 



(3vr),=o = l/g [|vr+)K°)|vr-) - |vrO)K+)|vr-) - + 

|7r-)|7r+)K°) - |7r-)|7r0)|7r+) + K0)|7r-)K+)] . (2.4) 



2.1. Decay of the 77 meson into vr tt tt 
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Figure 2.1: Kinematical boundaries of the Dalitz plot for the rj — > tt+tt 7r° decay. 

This wave function given above is antisymmetric under any exchange of pions: vr'^ -H- vr"*", 
Tr~ -H- 7r+ and ir^ -f-)- tt". In particular, by applying charge conjugation (equivalent to 
7r+ -f-T- TT^ exchange) one gets: 



This is in contradiction with Ac = +1 of the rj meson. Therefore, the decay r] — )• 7r+7r~7r'^ 
in the isospin state / = violates C symmetry. One can show that in this case also G 
symmetry is broken. The G operator is given by: 



where C stands for charge parity, and I2 denote the operator of the second component of 
the isospin. The eigenvalue of this operator is given by Xq = {—l)^Xc, thus = — 1 
for pions and Ac = +1 for the rj mesons. Hence, the rj — )• tt~^-k~7:^ in the / = does 
not conserve G symmetry. Therefore, the decay r] — t- vr+vr^vr'^ may occur if the isospin 
is conserved but then it has to violate C and G symmetry, or it may decay conserving C 
symmetry but then isospin is violated. 

As it was mentioned before this decay has a strong isospin breaking part which is 
driven by the term of the QCD Lagrangian proportional to the — m„. This isospin non- 
conserving interaction result in the final isospin state of three pions with 1=1 and Ac = 
+1. However, the interference between conserving and not-conserving charge symmetry 
production amplitudes of the Ac = +1 and Ac = — 1 states are also possible. This can 
result in asymmetry between kinetic energy distribution of charged pions vr"*" and vr". 




(2.5) 



G = Ce 



,inl2 
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Theory /Experiment 



A. 







d 



f 



34 



34 



ChPT LO 
ChPT NLO 
ChPT NNLO 
Gromley 
Layter | 26| 
Amsler 138 
Abele 
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34 



35 



36 



35 



120 
314 
538 
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Ambrosirii |i39| 



-1.039 

-1.371 

-1.271 
-1.18±0.02 
-1.080±0.014 
-0.94±0.15 
-1.22±0.07 



-1.090 



-0.013 
-0.024 



0.27 
0.452 
0.394 
0.20±0.03 
0.034±0.027 
0.11±0.27 
0.22±0.11 
0.124 ±0.016 



0.0 
0.053 
0.055 
0.04±0.04 
0.046±0.031 

0.06 (fixed) 



0.0 
0.027 
0.025 



0.057 



-0.013 
-0.022 



0.14 ±0.03 



Table 2.1: Dalitz plot parameters obtained from the theoretical predictions of the ChPT (first 
three rows), and the same parameters obtained from various experimental measurements of the 
77 — > Tr+Tr^TT*^ decay. 



A convenient way to study the ry — )• -k^tt^tt^ decay in view of the possible C violation 
is to use the Dalitz plot. For that purpose one can use the Mandelstam variables defined 
as: 



{Pv-Pif = ( 



TJlr 



2 • rrirjTi, 



(2.7) 



where pi and rm denote the four-momentum vectors and masses of final state particles, 
and Tj stands for the kinetic energy in the rest frame of the rj meson. However, in case of 
the 7r"''7r~7r'' final state where = , one can use the symmetrized and dimensionless 
variables defined as: 



X 



Y 



1, 



(2.8) 
(2.9) 



V Q 

31o 

Q 

where Q = T+ ± T_ ± Tq = — 2m^± — m^o is the excess energy. The kinematical 
boundary of the Dalitz plot in the X,Y plane for the r] — t- tt^tt^tt^ decay is shown in 
Fig.[2T} 

The distribution inside the boundaries is symmetric and fiat when the transition matrix 
element is constant. In general the density distribution is given by the matrix element 
squared which can be described by expanding the amplitude in the powers of X and Y: 



|M|2 = Al{l + aY + bV^ + cX + dX^ + fV^ + ...) 



(2.10) 



where a, b, c,d, f, ... are the parameters which can be obtained phenomenologically or on the 
ground of theory, and Aq stands for the normalization factor. By extracting the parameters 
from the experimental data and comparing them to the theoretical predictions one can test 
the assumptions of the Chiral Perturbation Theory. Furthermore, c coefficient, and other 
parameters standing in the odd-powers of X are sensitive to charge conjugation violation. 
The values of the parameters obtained from previous measurements and from calculations 




calculated based on the ChPT are collected in Tab. 2.1 In all previous experiments 



c parameter was found to be consistent with zero, but was not explicitly given in the 



publications due to large errors. Furthermore, only KLOE 39 experiment collected enough 
statistics to establish value of the f parameter. 

The amplitude mixing between Ac = — 1 and Ac = +1, describing the transition into 
isospin state 1=1 and / = 0,2, respectively, can be investigated by studying of the 
symmetries of population in different parts of the Dalitz plot. In particular the possible 
presence of C violation could be observed in three parameters: (i) left-right asymmetry 
- Alr, (ii) quadrant asymmetry - Aq, and (iii) sextant asymmetry - As- Each of these 
parameters depends on different isospin states of the final three pions. The asymmetries 
are defined as number of events observed in different sectors of the Dalitz plot divided as 



it is shown in Fig. 2.2 



The left-right asymmetry is defined as: 

A,n=''^-''\ (2.11) 

where the Nl stands for the number of events where vr^ has a larger energy than vr"*" and 
and denotes the number of events where the vr"*" has greater energy than vr^. It is 
sensitive to C violation averaged over all isospin states. However, it is possible to test the 
charge conjugation invariance in given I state. For this one uses the quadrant and sextant 
asymmetries which are defined as: 

_ N, + N,-N,-N, 

- N, + N, + N, + N,' ^^-^^^ 

^ Ni+N2 + N-i + Ni + N5 + NQ' ^ ■ ^ 

where Ni denotes the number of observed events in z-th sector of the Dalitz plot. The 
quadrant asymmetry tests the C invariance in transition into the Stt final state with / = 
2, and the sextant asymmetry is sensitive to the 1=1 |40]. Table 2.2 summarizes all 
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experimentally measured values of three asymmetries. The previously measured values of 
the asymmetry parameters indicate no C violation with respect to calculated uncertainties. 



Experiment 


Alr X 10-2 


Aq X 10-2 


As X 10-2 


Layter \26 
Jane 41 
Ambrosino 
PDG average 


39 
12 


-0.05 ±0.22 
0.28 ±0.26 
0.09 ±0.10 
-0 09+°-^^ 


-0.07 ±0.22 
-0.30 ±0.25 
-0.05 ±0.10 
-0.09 ±0.09 


0.10 ±0.22 

0.20 ±0.25 

0.08 ±0.10 

l2+°-i° 
'-'•J^^-o.ii 



Table 2.2: The values of Dalitz plot asymmetries obtained experimentally together with the 
average value from the PDG for the rj Tr+vr-Tr" decay. 



2.2 Decay of the rj meson into n^e^e 

The investigation of the charge conjugation invariance C in the electromagnetic inter- 
actions can be done by studying the r] — t- vr'^e^e- decay. In the framework of the Standard 
Model and the QED the matrix element for this process should involve the two virtual 
photon exchange [42^ as it is presented in Fig. 2.3 with the transition according to the 
reaction: 

r/^7r°±7*±7* ^7r°±e+±e-. (2.14) 

Therefore, the wave function of the 7r'^7*7* system transforms with the C operator as 
follows: 

C(A*7*) = A^,"aJ AJ'(A*7*) = (+1) • (-1) • (-1)(A*7*) = ±1(A*7*). (2.15) 
The eigenvalue of the charge parity for the rj meson is Ac = ±1 which is in agreement with 



^ > 




Figure 2.3: Decay mode r] — > 7r"7*7* — s- 7r°e+e occurring by the C-conserving second order 
electromagnetic process. 

above shown charge parity of the decay system, thus this process hold the C invariance. 



2.2. Decay of the rj meson into iv'^e'^e 
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Figure 2.4: Diagram for the C invariance violating transition 77 — > n'^'y* — >■ Tr^e+e occurring by 
the first order electromagnetic process. 

The decay rate of this C-conserving process, predicted theoretically ranges from 10^^^ to 
10~^ depending on the undertaken assumptions: 

BR{r] 7r°e+e~) w (1.5 ± 0.4) • 10"^^ 
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BR{r] 7r°e+e") « 1.1 • 10 



-8 
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BR{ri ^ TT^e+e-) ^ (1-6) -IQ-^ [45]. 



(2.16) 

(2.17) 
(2.18) 

It is worth to mention that the second and third predictions are based on the approach of 
the Vector Dominance Model (VMD) [46^. In the framework of this model it is assumed 
that the decay is dominated by the virtual transition r] — )• Vj* followed by the V — )• 
and 7*7* — )• e~^e~ , where the V denotes all the neutral vector mesons of zero strangeness: 

However, in principle the decay t] — )• 7r^e'*"e~ may also be realized in first order of 
the electromagnetic interaction with only single 7 quantum in the intermediate state (see 



Fig. 2.4) via transition: 



7? 7r° + 7* 7r° + e+ + e . 



(2.19) 



In this case the C operator acting on the wave function of the intermediate state leads to: 

C(A*) = AJ.''Aj*(vrV) = (+1) • (-l)(vrV) = -1(A*), (2-20) 

which is in contradiction with the charge parity of the rj meson (A^ = +1). Thus, this 
process introduces the violation of charge conjugation. The decay width for the first 
order electromagnetic processes are larger than for second order mechanism. Therefore, 
experimentally the C-invariance breaking would then manifest itself with increasing the 
branching ratio with respect to the predictions listed in equations 2.16, 2.17 and 2.18. 
At present only an experimental upper limit for the rate of the branching ratio BR{r] — )• 

, and it amounts to 4 • 10~^ fl2|. Therefore, still at least 



was determined 
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three orders of magnitude remains to be experimentally investigated until value predicted 
based on the Standard Model will be reached. The observation of higher branching ratio 
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than one calculated in the framework of the Standard Model could provide the evidence 
that the decay r] — )• Tr^e~^e~ is not conserving C-invariance. 

One has to stress that C- violating first order electromagnetic process is the most prob- 
able possibility, but not the only one. Another conceivable process which may lead to the 
increase of the branching ratio could be e.g. rj — )• x^"^^ ~^ j*^*7r^, where is an unknown 
particle not included in the Standard Model |53| . 

The aim of this thesis is to contribute in searches of the C violation by determining the 
BR[r] — >■ Tr^e~^e~) or lowering the present upper limit for this branching ratio. 



3. Experimental methods 



The measurement described in this thesis was carried out in October and November 
2008 by means of the Wide Angle Shower Apparatus (WASA) |54| and the Cooler Syn- 
chrotron (COSY) [55] operating at the Research Center Jiilich, Germany. The rj meson 
was created in proton-proton collisions via pp — t- pprj reaction with the beam momentum 
oipbeam — 2.142 GeV/c which corresponds to the kinetic beam energy of Tf,eam = 1.4 GeV. 
The experiment was based on measurement of four-momentum vectors of outgoing nucle- 
ons and of decay products of unregistered short lived t] meson which was identified using 
the missing and invariant mass techniques. 



3.1 Cooler Synchrotron COSY 

The COoler SYnchrotron "COSY" is a storage ring which can deliver unpolarized 
and polarized proton and deuteron beams in momentum range between 300 MeV/c and 
3700 MeV/c. The ring consists of 24 dipoles and 56 quadrapole magnets which are used 
to keep and focus particle trajectories during the acceleration process, and also sextupole 
magnets which are used to deflect beam what results in achieving better beam optics. The 
acceleration process takes place in two steps: (i) first the ions {H~ or D~) are accelerated 
in the isochronous cyclotron (JULIC) and next (ii) the beam is stripped of electrons and 



finally injected into the 184 m long COSY ring (see Fig. 3.1), where particles are stored 
and accelerated up to the demanded momentum. The beam is then directed into internal 
or external experimental targets. The beam energy range allows for production of all basic 
pseudoscalar and vector mesons up to mass of 0(1019) particle. 

Additionally COSY accelerator is equipped with two types of beam cooling systems: 
(i) an electron and (ii) stochastic, used for low and high energies, respectively [56]. Both 
cooling methods allow to reduce the momentum and spatial spread of the beam. The 
beam momentum spread ^ after applying both types of cooling method can be reduced 
to around 10~^ 57 . 



The whole process of acceleration together with the beam cooling phase takes a few 
seconds. COSY storage ring can be filled with up to 10^^ particles, and the life time of 
the circulating beam varies from minutes to hours depending on the thickness of used 



target. Presently in COSY accelerator facility two internal experiments: WASA 54 and 



ANKE \58] are in operation, and one external experiment: COSY-TOF 59 
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Figure 3.1: Schematic view of the COSY Cooler Synchrotron 55 storage ring at the Research 
Centre JiiUch. Presently used detector systems: WASA |54l, COSY-TOF (on external beam 
line) (59] and ANKE (58], are marked in yellow. 



3.2. WASA-at-COSY apparatus 
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The Wide Angle Shower Apparatus - WASA detector - originally operating at the 



CELSIUS 60 61| facility in Uppsala, Sweden was transferred to COSY accelerator facility 

is a 



in 2006 62 . The new WASA-at-COSY 54 detector, shown schematically in Fig 3.2 



large acceptance detector consisting of three main parts: the Central Detector (CD), the 
Forward Detector (FD) and the Pellet Target system. 



Pellet line 



TOF Detector 



Solenoid 




EM Calorimeter 



Plastic Scintillators 

Range Hodoscope 



Figure 3.2: Schematic cross view of the WASA-at-COSY apparatus. The detector components 
are described in the text. 



WASA-at-COSY detector system is capable to register neutral and charged particles 
emerged in the collision of beam and target and also particles originating from the decays 
of short lived mesons. It was developed mainly for production and detection of the vr'^ and 
rj meson decay products in order to study the fundamental symmetries and to test the 
Standard Model. 



3.2.1 Central Detector (CD) 

The WASA Central Detector (CD) |63[[64] is positioned around the beam and target 
interaction point. It is used to detect and identify light neutral and charged particles like: 
7, e"*", e~ , TT"^, 7r~ which originate from decay of the short lived mesons and direct collision of 
nucleons. The most inner part of CD is the Mini Drift Chamber (MDC) surrounded by the 
Plastic Scintillator Barrel (PSB) and the yoke of the super-conducting solenoid |65| which 
together enables to measure momenta of charge particles. The outer part constitutes the 
Scintillating Electromagnetic Calorimeter (SEC) which is used to measure particles energy. 
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Mini Drift Chamber MDC 

The MDC detector is mounted around the beam pipe inside the superconducting 
solenoid which provides an axial magnetic field up to B = 1 T. It is arranged in a cylindrical 
shape consisting of 1738 straws stacked in 17 layers (see Fig 3.3). Each straw tube is made 
out of thin (24//m) mylar foil, aluminized from inside, and with a gold plated sense wire 
in its center with diameter of 20fim 66 . The MDC has nine layers with straws parallel 
to the beam direction and eight with a small skew angle with respect to the beam line. 
The tubes are filled with gas mixture: 80% of argon and 20% of ethane to ensure that 
each charged particle which pass through a single tube will cause an ionization. The MDC 
vertex detector enables to determine the momenta of charged particles |67| with accuracy 
of Ap/p < 1% for electrons and positrons, Ap/p < 4% for charged pions, and Ap/p < 5% 
for protons i68j. 



Plastic Scintillator Barrel — PSB 

The PSB detector is used to determine the energy loss of charged particles and particle 
identification by AE — |pj method. The 52 scintillator bars of PSB detector are arranged 



in cylindrical shape 69 70] around the straw drift chamber where each bar is overlapping 



with next one with 7° to assure the coverage of a full geometrical acceptance (see Fig 3.3 



(right)). Additionally the forward and backward part is equipped with "end-cups" made 
out of trapezoidal elements arranged around the beam pipe. 



Scintillating Electromagnetic Calorimeter — SEC 

The Scintillating Electromagnetic Calorimeter consists of 1012 sodium doped cesium 
iodide CsI(Na) scintillating crystals in a shape of the pyramids in order to be arranged 



3.2. WASA-at-COSY apparatus 
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size (cm) 




Figure 3.4: (left) Schematic 3D view of the Scintillating Electromagnetic Calorimeter (SEC) 
consisting of 1024 scintillating modules arranged in 24 layers, (right) Angular coverage of the 
Scintillating Electromagnetic Calorimeter (SEC) . 



spherically around the interaction point (see Fig. 3.4 (left)). The CsI(Na) scintillating 
material provides a large light yield and short radiation length making them a very good 
material for measuring the energy and scattering angles of neutral and charged particles 
such as gamma quanta, electrons, positrons and pions 64 , 7l]. The crystals are grouped in 
24 layers covering almost the 47r acceptance (in polar angle: 20° < 9 < 169°, and azimuthal 
angle: 0° < (j) < 360°). Depending on the layer the length of the crystals varies from 30 cm 
in central part, to 20 cm in backward and 25 cm in forward part (see Fig. 3.4 (right)). 
The overall energy resolution of the SEC for photons can be described by the relation: 
^^^^ = ^= and the angular resolution for scattering angle is equal to 5° (FWHM) 



64 



The "punch-through" kinetic energy for pions is 190 MeV and for protons 400 MeV. The 
electrons, positrons and photons are stopped in the calorimeter depositing all its energy. 



3.2.2 Forward Detector (FD) 

The Forward Detector (FD) of the WASA apparatus consists of fourteen scintillating 
layers and a proportional straw drift chamber, and its geometrical acceptance covers in 
laboratory frame a range in polar angle from 3° to 18°. Such a setup enables to measure 
the energy loss and trajectories of recoil particles, mainly protons, deuterons and ^He 
nuclei. The particle identification in FD is based on the AE — E method which enables 
to reconstruct proton energy with overall resolution of about 10%. To further improve 
energy resolution and particle identification a new reconstruction technique based on the 



Time- Of- Flight measurement is under development [15,72-75 . Another possible upgrade 
which is in progress is a Cherenkov DIRC detector |76[]78| . 
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Figure 3.5: (left) Forward Window Counter (FWC) build out of two layers each consisting of 
24 cake-piece modules, (middle and right) Forward Proportional Chamber (FPC) build out of 
4 layers of straw tubes. 



Forward Window Counter — FWC 

The Forward Window Counter (FWC) is a closest detector to the conical exit window 
of the axially symmetric scattering chamber. The detector is 48-fold segmented and it 
is composed of two layers a 24 cake-pieced elements made out of 3 mm thick plastic 



scintillator 179 80 . The first layer is arranged in a conical shape whereas the elements of 



the second layer are assembled in a vertical plane (see Fig. 3.5 (left)). The elements of 
the second plane are rotated by one half of module - 7.5° - with respect to the first layer. 
Such an arrangement ensures a complete coverage of the forward scattering area. The light 
collection in the hodoscope is optimized to keep the detection efficiency as homogeneous 
as possible over the full detector. It is worth to mention that this detector will serve as a 
"start" detector for the Time-of-Flight method. 



Forward Proportional Chamber — FPC 

The Forward Proportional Chamber (FPC) is a detector located directly after the FWC 
and it is used for trajectory reconstruction purpose. It can measure particles scattering 



angles with a precision better than 0.2° 81,82 . It consists of 1952 thin straws stacked 



by 122 in sixteen layers and grouped in four detection modules. Each module is turned 
with respect to each other by 45° |83|. With respect to the x-axis the angles position of 



subsequent layers are: 315, 45, and 90 degrees (see Fig. 3.5 (right)). Each straw has a 
diameter of 8 mm and it is made out of mylar foil aluminized from inside with a 20 fim 
stainless steel sense wire in its center. All straws are filled with a gas mixture: 80 % of 
argon and 20 % of ethane to ensure an efficient ionization. 
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Forward Trigger Hodoscope — FTH 

The Forward Trigger Hodoscope (FTH) consists out of 96 individual plastic scintillator 
modules arranged in three layers: (i) two layers with a 24 elements each, in a form of 
Archimedean spiral rotated clockwise and counterclockwise, (ii) and one layer with 48 
cake-piece shaped elements |84[j86| . Overlap of these three layers gives 24 x 24 x 24 pixel 



map (see Fig 3.6). Whole detector setup of the FTH has highly homogeneous detection 
efficiency and shows a fairly uniform behavior jsTj. The FTH is used as a first level 
trigger and gives information about particles multiplicities. It is also possible to use FTH 
to real time scattering angle reconstruction of individual tracks and combine it with the 
information about deposited energy in successive layers and to use it to determine the 
missing mass of forward going particles on line on the trigger level (87|. FTH can be also 
used to determine the energy losses and thus can be used as AE detector for identification 
of recoil particles via AE — E method. 



Forward Range Hodoscope — FRH 

The Forward Range Hodoscope (FRH) is build out of five thick plastic scintillator 
layers, each cut into 24 cake-piece elements [88|. The first three layers have thickness of 



11 cm while layers 4 and 5 have 15 cm (see Fig. 3.7 (left)) [89] . The FRH enables to 
reconstruct kinetic energy of charged particles and to use the AE — E method for particle 
identification. The relative kinetic energy resolution for protons up to 360 MeV which are 
stopped in FRH, is almost constant and is equal to about 3 %. For more energetic protons 
resolution worsens linearly to be about 10 % for protons with energies about 600 MeV. 
FRH provides also information for the trigger matching algorithm to verify alignment of 
each track in the azimuthal plane. 
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Figure 3.7: (left) View of the Forward Range Hodoscope (FRH) build out of the cake-pieced 
scintillating modules arranged in five layers, (right) Front view of the Forward Veto Hodoscope 
(FVH) with 12 horizontal scintillation modules readed from both sides. 



Forward Veto Hodoscope — FVH 

The Forward Veto Hodoscope (FVH) consists of 12 horizontally and 12 vertically placed 



plastic scintillator bars equipped with photomultipliers on both sides 80 90] (see Fig. 3.7 



(right)). This enables to reconstruct particle hit position from time signals registered on 
both sides of the module. The bars are arranged in two layers with the relative distance 
of 77 cm. The FVH is mainly used to detect particles punching through the FRH and 
to reject them as too energetic for the pp — t- ppij reaction. Depending on the measured 
reaction, a passive iron absorber can be placed between FRH and FVH. The thickness of 
the absorber can be chosen from 5 mm up to 100 mm. Usage of the absorber enables to 
disentangle between slow and fast recoil protons coming from meson production reaction 
and elastic scattering. 



3.2.3 Pellet Target system 

The WASA-at-COSY detector is equipped with a specially designed target system [9l| - 
[94] , providing high density frozen droplets of hydrogen or deuterium, called pellets. It is 
located on a platform above the Central Detector, delivering pellets of an average size of 
35 /im, with frequency rate of about 10 kHz, to the interaction region by a thin 2 m long 
pipe. The droplets are produced by a piezoelectric transducer, which induce vibration of 
the nozzle and brakes liquid stream into pieces. After passing through the scattering cham- 



ber, pellets are captured in a cryogenic dump. Fig. 3.8 shows schematically the target-beam 
arrangement (left) and the scheme of the pellet generation process. Such a construction 
satisfies requirements for achieving high densities of particles up to 10^^ atoms/cm^ re- 
sulting in luminosities up to lO^'^ cm~^ s~^ with combination of the 10^^ particles of the 
COSY beam. Thus, WASA-at-COSY is capable to carry out high statistics experiments, 
needed to study rare and very rare decays of mesons. 



3.3. Data Acquisition System and trigger logic 
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Figure 3.8: View of the WASA-at-COSY Pellet Target System 91-94 



3.3 Data Acquisition System and trigger logic 

In WASA-at-COSY experiment, to handle a high event rates which are typically larger 
than 10 kHz, an efficient Data Acquisition System (DAQ) was developed to digitalized and 
store data for further off-line analysis [95] . The schematic view of the system is show in 
Fig|3;9] 

Variety of detectors used for registering desired processes, forces to use different dig- 
italization modules. The plastic scintillator elements (FWC, FTH, FRH, FVH and PS) 
are read out by photomultipliers, from where the signal is directed to the splitters which 
divide it (i) and directs to the leading edge discriminators and then to the trigger system, 
and (ii) to the shaper which stretches the signal to about 100 ns and then to the Fast- 
Charge-to-Digital-Converters (F-QDC) which converts analog signals to digits (numbers) 
which can be further processed by computers |96| . Also the same signal is used to digi- 
tize time information by Fast-Time-to-Digital-Converters (F-TDC-GPX) which has a time 



resolution of about 85 ps 97 



The signals from the Electromagnetic Calorimeter crystals are splitted into two branches. 
One goes to readout system built out of Slow-QDC's (Flash-ADC chips) and then to the 
Field-Programmable-Gate-Array (FPGA) to be integrated for charge measurement. The 
Slow-QDC's can run in "floating gate" and "fixed gate" mode depending on the experimen- 
tal demands. Second signal is sent to the discriminators and used for summing different 
groups of signals together and then to be applied as a logic signals which are sent directly 
to the trigger system. 

The straw chambers (FPC and MDC) signals are ffist amplified and then sent to the 
discriminators and then they are digitalized by the Fl-TDC's (Slow-TDC), with the time 
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Figure 3.9: Block scheme of the Data Acquisition System (DAQ) of the WASA-at-COSY exper- 
iment. 



resolution of around 120 ps 

The whole system consists of 14 crates equipped with QDC's and TDC's digitalization 
modules, which continuously samples data streams and stores signals from each module in 
FIFCQ queue. This type of buffer allows to run the measurement without a trigger 
delay. Each one of the digitalized signal is marked with a time stamp relative to the trigger, 
which is broadcasted by a Synchronization System connected to each crate by Low- Voltage- 
Differential-Signaling (LVDS) bus. After the digitalization, signals are matched by the time 
stamp and marked with the same event number. 

Before event can be saved it has to be checked to fulfill the trigger logic j99j. The 
conditions of the trigger are based on the time and geometrical coincidences and hit mul- 
tiplicities in different detection modules. In the Forward Detector it is also possible to 
apply a simple track finding algorithm on a trigger level. The so called "Matching Trigger", 
compares hit position in consecutive layers of FWC, FTH, FRH and decides if they are 
coming from the same particle 100 . This technique allows to select interesting events dur- 
ing the measurement and reduce the data rate to be later stored on the disk. Depending 
on the studied physical process several trigger conditions can be imposed simultaneously 
on the event by applying logical "and" operation. Also a multiple trigger conditions can 
be applied to the same data stream at the same time. 

Further on, after passing through the trigger level, events are sent to Event Builder 



^An acronym for First In, First Out, an abstraction related to ways of organizing and manipulation of 
data relative to time and prioritization. 
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which stores them in files of 20 GB size on a RAIE|^ system. Each file is marked by a 
unique run number and after some time it is transferred from RAID to the tape archive 
for permanent storage. 



3.4 Missing and invariant mass techniques 

In order to evaluate physical observable like branching ratios and asymmetry param- 
eters, one has to extract clean signal of desired reaction from the measured data sample. 
In this thesis two complete reaction chains have to be identified: 

pp — )• ppr] — )• ppn^TT^ir^ — )• ppn^TT^jj, (3.1) 

pp — )• ppr] — )• pp7r°e"'"e~ — )• ppe^e~'~f'~f. (3.2) 

In both cases reaction identification will relay on determining the four-momentum vectors 
of all particles in final state. The reconstruction of protons in FRH will be based on 
measurement of energy loss and the direction vector r in the FPC. The charged pions 
and electrons will be registered in CD, where MDC will provide the information about 
the momentum vector p, and energy losses will be measured in PS and SEC. The gamma 
quanta originating from the decay of neutral pion will be registered in SEC and their, 
four-momentum vectors will be reconstructed based on the energy losses and the positions 
of the hits. 

The procedure of extraction of signal from interesting reaction is divided in several 
analysis steps, where conditions for the minimal thresholds of energy are set and time 
coincidences between different tracks are checked. In general the signal identification will 
be based on the reconstruction of missing and invariant masses. The missing mass is 
defined as follows: 



{El - p-^) = J {Eh + mt- Ep, - Ep,f - {n - pi - p2)2, (3.3) 



where the rux denotes the mass of unregistered short lived meson (in our case rj meson), 
Sfe, ph corresponds to the energy and momentum vector of the beam, respectively, and Ep^ , 
^P2i pi) P'i represent energies and momenta of recoil protons, and the mt stands for the 
mass of the target. While the invariant mass reads: 




where Ei and pi corresponds to the energies and momenta of decay products of short lived 
meson. 



■^An acronym for Redundant Array of Independent Disks which is a storage technology that combines 
multiple disk drive components into a logical unit. 
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As a first stage of reaction identification a missing mass of tlie system pf,pt — )• 

Next cliarge particles (7r^,e^) will be 



will be determined according to the equation 
identified via invariant mass: 



3.3 



7r±(e±) Pn±{e±)- 



(3.5) 



Finally the signal from the vr*^ meson will be reconstructed as the invariant mass of two 



gamma quanta according to the equation 3.4 



4. Reaction kinematics, analysis 
methods and detector simulation 



In order to understand properly the physical processes, the computer simulations of 
the detector response for both investigated reactions were carried out. For that purpose 
we used the hadronic event generator Pluto++ [lOl] and a WASA Monte Carlo Software 
which is based on the GEANT package 102| used to generate the response of whole detector 
setup. Based on this we were able to reproduce the kinematics of the investigated reactions 
and to calculate the geometrical acceptance of the detector. 



4.1 Event generator: PlutoH — h 

The simulations were performed using the Pluto++ event generator, which enables 
to simulate kinematics of reactions at beam momentum given by the user. As an output 
the four-momentum vectors of all particles in the final state are returned. In most of the 
cases, four-momentum vectors of all particles in the exit channel are weighted according to 
the Phase Spac^ using GENBOD routines 103 . For a few reactions and decay channels 



a phenomenological properties of a given process like angular distributions and transition 
matrix elements are implemented into the code. 



4.2 WASA Monte Carlo package 

To simulate the response of the detector for the generated data, a virtual model of 



the WASA detector was built based on a GEANT 3.21 framework 102 . The description 
of the detector elements and support structures were fully modelled using an abstract 
objects called "volumes" which are filled with appropriate materials and embedded in a 
3D coordinate system with an interaction point in its center. Each volume is described by 
the geometrical dimensions and positions relative to the center of the coordinate system. 
The implementation includes the sensitive, as well the passive materials which creates the 
whole apparatus. 

The interactions of the particles with the detector is done by propagating each of them 
through the volumes and simulating according to the known cross sections the physical 
processes like: photon conversion, production of secondary particles, quenching effects. 



^In this context the Phase Space term refers to an isotropically and uniformly distributed values of the 
four-momentum vectors of particles produced in the 's' wave with the relative angular momentum J=0. 
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Figure 4.1: The schematic diagram of the Root-Sorter analysis workflow. 



multiple Coulomb scattering, hadronic interactions, energy losses. However, the package 
do not include the light propagation processes in the scintillator, electronic noise, and the 
response of the photomultipliers. But still, the package provides a very precise comparison 
to the measured experimental data. 

The Monte Carlo simulation starts with the input of four-momentum vectors from 
Pluto++ event generator, and is continued with the stepping action of each particle through 
the volumes. Generated responses for each event include the map of activated detector 
elements, values of deposited energy and times. Further, matching of actual experimental 
conditions with the simulations, an additional smearing of the observable is possible via 
separate input filter files. 



4.3 Root-Sorter work flow scheme 

The main analysis software used for off-line data processing is the Root-Sorter frame- 



work [104j which bases on the object oriented C++ language [105 and CERN-ROOT 



package 106 . The software is organized in a modular way, such that, each class is re- 
sponsible for different tasks like: tracks or cluster finding, calibration, and reconstruction. 
The ROOT software delivers very rich library of standard mathematical functions, his- 
togramming and fitting procedures, and also provides the management of the input/output 
streams. It provides very useful classes for the vector and four-momentum vectors ma- 
nipulations as well most of the non-standard and dedicated methods of reconstruction. 
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calibration and track finding which were prepared by the members of the WASA-at-COSY 
and CELSIUS/WASA collaborations. 

The event processing is divided into two stages: (i) low level, and (ii) high level analysis 
(see Fig. 4.1). In the low level analysis, the raw signals from stored data files are loaded 
into the computer memory and then calibrated from QDC and TDC units, to the units 
of energy represented in gigaelectronvolts (GeV) and time expressed in nanoseconds (ns). 
Also, the offset due to the electronic delays and non-uniformity corrections are applied on 
this stage. For the Monte Carlo data, instead of the calibration procedures, the filters are 
applied, which stores the smearing values of energy and times in order to account for the 
experimental resolution of the detection system. 

The data flow for experiment after decoding and calibration, and for simulations after 
filtering is organized in the same way. At the beginning the signals from individual detec- 
tion modules are composed into hits, and then they are stored in a HitBanks separately 
allocated for each detector. Then the cluster algorithm groups the adjacent hits coming 
from one detector into clusters which are stored in Cluster-Banks. Subsequently, the track 
finding algorithm is applied to the cluster banks and creates from the clusters belonging 
to different detectors an object, which is representing a particle track. 

The high level analysis refers in most of the cases to the stage, in which user by himself 
writes an analysis procedures to extract signal of desired reaction, using reconstructed 
tracks and constructing missing an invariant mass spectra, and also applying cuts. On 
this stage user loads the low level analysis modules and, the event processing is then 
automatically performed, so that the user can access information about registered particles 
in the form of the deposited energy and track direction r. Furthermore, for backward 
consistency, the track-, cluster- and hit- banks have an ability to inherit properties from 
its predecessors, so that it is possible for the user at any stage of the analysis to access to 
low level information from top-level data structures. 



4.4 Kinematics of the pp — > ppr] — > ppn^n n^i'yj) 
and pp — )■ ppr} — )► pp7T^{'y^)e^e~ reactions 

In the experiment the production of the r/ meson was performed by colliding the 
proton beam of 1.4 GeV energy with a stationary proton pellet target. For the simulation 
the Pluto++ event generator was used which includes properties of the rj meson production 
mechanism. One of them is the anisotropy parameter of the rj polar angle in the center of 
mass frame 0^"^. This anisotropy was observed by the DISTO Collaboration 



107 



at the 



beam kinetic energies of 2.15, 2.5 and 2.85 GeV. In case of the investigated reaction the 
beam energy corresponds to the excess energy of Q = 55 MeV. As it was shown by the 



COSY-11 and COSY-TOF collaborations 109-112 , in reactions closer to the kinematic 



threshold for the r] meson production this anisotropy vanishes. In addition, the proton 
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Figure 4.2: The simulated kinematics plots shows the correlation of kinetic energy with the 
scattering angles for particles in the final state: (upper left) protons, (upper right) gamma 
quanta, (lower left) pions, (lower right) electrons and positrons. The red lines indicate the 
geometrical acceptance of the WASA-at-COSY detector. 



angular distribution, which with the increasing beam momentum shows tendency to be 
aligned stronger to forward and backward directions was included. The decay products of 
the 77 meson were distributed homogeneously in the whole phase space. For both reactions 
the distribution of the angles and energy of recoil protons and gamma quanta originating 
from the tt^ decay, are the same. The differences appears in the spectra of the charged 
particles. As it is shown in Fig. 4.2 (lower left) the pions vr^ populate more frequently 



the lower kinetic energy region and lower polar angles. While, the electrons and positrons 



have more broader distribution of the energy and 9 (Fig. 4.2 (lower right)). This difference 
is due to around 300 times lower mass of the electrons than pions. 

In order to estimate the acceptance, the geometrical coverage of the CD and FD was 
taken into account. For detecting protons in forward direction the acceptance is equivalent 
to the angular range of 3° < 9 < 18° and for pions and gamma quanta in central part is 



4.4. Kinematics of the pp — >■ ppi] — >■ ppiv'^TT Tv'^ij'f) and pp 



PPV 



pp7r''(77)e"'"e reactions 41 



20° < < 169° (see Fig. 4.2). For the reaction pp — )• pprj — )• vr+vr 'K^{'yj) the fraction of 



events which were found in the sensitive range of the detector (red hnes) is equal to 35.3%. 
While, for the pp — )• ppr] — )• 7r''(77)e'*"e~ process the geometrical acceptance is equal to 
46.9%. 

To check the event generator properties the Dalitz plot expressed in variables defined 



by equations (2.8) and (2.9) was plotted, using the simulated four-momentum vectors. As 
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Figure 4.3: Simulations of the Dalitz plot for the ?/ tt+tt^tt^ decay according the to the 
homogeneous phase space population, (left) The two dimensional Dalitz plot in normalized X,Y 
coordinate system, (middle) Projection of the Dalitz plot to X axis, (right) Projection of the 
Dalitz Plot to Y axis. 

it was mentioned in Chapter 2 the Dalitz plot represents a population of kinematically 
available phase space. In the case of simulation, performed under the assumption of no 
interactions between particles, the Dalitz plot should be populated homogeneously. As it is 
shown in Fig. |4.3| (left), the results of the simulations are consistent with these expectations. 



Additionally in Fig. |4.3| (middle) and Fig. 4.3 (right) the projections of the Dalitz plot into 
the X and Y axes are presented. 
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5. Track selection and reconstruction 



The data analysis aiming at identification of pp — )• ppr] — t- pp7r+7r^7r*^(77) and pp — )■ 
pp-q — )• pp'rT^{'^'^)e^ e~ reaction cliains, at tlie first stage is tlie same, tlius tlie procedures 
for botli reactions, will be discussed simultaneously. The track reconstruction is based on 
signals registered in Forward and Central Detector of WASA setup. At the beginning of 
the reconstruction process hits in detector elements which belong to the same particle are 
merged into clusters. Next, the clusters reconstructed in different detectors are combined 
into single particle track. For both investigated processes the trigger conditions and track 
reconstruction procedures are the same. 



5.1 Data set and trigger conditions 

The data analyzed in this thesis were collected during the two week run in the year 
2008. The r] meson was produced in a reaction of proton beam with momentum of pi, = 
2.142 GeV/c and a stationary proton target. The excess energy for the ppr] system was 
equal to Q = 60 MeV, for which the production cross section of the r] meson amounts to cr = 



9.8± 1.0 nh 113). This cross section allowed to measure large event rates enabling to study 



rare and very rare decay processes. However, in the proton-proton interactions one has 
to deal with a large background originating from the direct multi-meson production. For 
the interesting 77 — t- vr^Tr^vr'^ decay the physical background comes from direct production 
of three pions via the pp — )• ppir^-K^TT^ reaction channel. The total cross section for 
this reaction is in the same order of magnitude as the one for the 7] production process. 
Therefore, part of this background with an invariant mass of Tr+vr^vr'^ system close to the 
mass of the rj meson will contribute to the signal range. However, expected signal to 
background ratio should be about ten in case of tagging by means of the missing mass 
of two protons, with resolution of a few MeV |114 . Yet, the situation for a two pions 



direct production {pp — t- ppir^iT") is worser, because the total cross section is hundred 
times larger |115| than for the r] meson production. But in this case only events with 
misidentified Tr*^ can be mistakenly taken as the signal and therefore one can suppress the 
background due to different final state than signal reaction. 

In case of the rj — t- ■iT^e~^e~ final state, the background comes mostly from the reaction 
rj — )• e'^e^'j where the additional gamma quantum is reconstructed due to the splitting of 
signals in the calorimeter, and reactions with final state of tt^tt^tt^ for which the electrons 
and positrons can be misidentified with charged pions. Another source of background 
constitutes the conversion of gamma quanta on the beam pipe, where the interaction 
of photons in the beryllium can cause emission of the electron and positron pair. To 
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Figure 5.1: The ADC distribution for the Forward Veto Hodoscope detector after pedestal sub- 
traction. The threshold on the discriminator used for the T^Yq^ was set to 22 mV which corresponds 
effectively to the value of 300 in the units of the ADC channels. 

suppress this type of background the studies of the conversion process has to be performed. 
Moreover, the direct production of two neutral pion with subsequent decays, into vr*^ — )• 77 
and vr'^ — )• e~^e~^ {pp — )• vr'^vr'^ — )• e'^e~777) can obscure the searched signal. Thus the 
experimental conditions and large hadronic background make the investigation of the r] 
meson decays in proton-proton collisions experimentally challenging, and demands a very 
selective trigger. 

In order to reconstruct 2p, 27r (or 2e) and 2-f one has to construct a trigger, which 
basing on the simultaneous fulfilment of specific conditions in FTH, FRH, FVH, PSB and 
SEC, will select two charged particles in the Forward Detector and two charged and two 
neutral particles in the Central Detector. The first requirement was that at least two 
charged particles in second layer of the FRH will be detected T^^^^'^ (where jj. indicates 
the hit multiplicity in the detectors), together with the Matching Trigger condition T^i^ 
requiring that at least two groups of clusters from different detectors corresponds to the 
same azimuthal angle. Furthermore we set a rejection 'veto' condition for events when at 
least one charged particle was registered in FVH detector {T^^Yq^). This requirement was 
used to discard events with too energetic protons coming most probably from the elastic 
scattering reaction. The corresponding ADC spectrum of the FVH detector is shown in 



Fig. 5.1 In the Central Detector at least two charged particles in Plastic Scintillating 
Barrel T^^.^ were required and at least one particle in the Scintillating Electromagnetic 
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Calorimeter T^^f^. Thus summarizing, the main trigger for detection of the interesting 



reactions, reads: 

rp _ rpFRH2 . rpMT . rpFV H . rpPSB . rpSEC / c 1 ^ 

All the events which fulfilled this requirements were saved to the disk. 

5.2 Selection of tracks in the Forward Detector 

Finding and reconstruction of particle track in the Forward Detector relay on geo- 
metrical overlap of clusters, which are formed based on hits in different detectors. The 
place where a particle crosses the FTH, is determined as an overlapping region of three 
elements from each detection layer. This allows to determine a line from an assumed in- 
teraction point of the beam and target to the center of the pixel. Also a time information 
is determined for tracks as an average time of signals in each layer of FTH. Then, the 
angular information is improved by reconstruction of tracks in FPC, where the polar {9) 
and azimuthal {(f)) angles of the track are reconstructed based on position of sense wires 
and the drift times. For the reconstruction procedure a line parameters determined in FTH 
are used as an initial values. This procedure typically improves the angular resolution by 
a factor of two and more. After the track formation, the information from clusters about 




X[cm] 



Figure 5.2: (left) Schematic view of the FTH pixel formation in the XY plane perpendicular to 
the beam line. The green point denotes the center of the pixel which is taken as the hit position, 
(right) Example of the hit distribution in the FTH. The spectrum shows the hits in the center of 
each pixel. 

azimuthal overlapping, energy deposits, and time differences from the FWC, FRH and 
FVH are taken into account. 
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5.3 Selection of tracks in the Central Detector 

The Central Detector allows to register charged and neutral particles coming from the 
rj meson decays. As a first step, signals registered in each of the detection elements are 
formed into a cluster, for each detector separately. Next, clusters from different detectors 
are matched in order to reconstruct a track of particle which passes through the detection 
system. The track reconstruction algorithm, based on the signals registered in MDC and 
PSB searches for all charged particle tracks. Afterwords, remaining signals in SEC not 
associated with a charged particle, are treated as originating from neutral particles. 

In SEC when a particle hits the crystal an electromagnetic shower is produced, which 
spreads over nearest elements. The reconstruction procedure starts with a search for a 
crystal with deposited energy greater than 5 MeV. When such a crystal is found, the 
routine puts it in a center of a pre-cluster which is a square of 3 x 3 elements. Next, the 
routine checks energy in the adjacent elements. If the energy is greater than 2 MeV the 
element is added to the cluster. These elements also become a reference in the next step 
of the procedure. Finally, the algorithm terminates if there is no more crystals with the 
energy greater than 2 MeV. In addition only elements which have a time signals within 
50 ns window are included into a cluster. Further on, only clusters with sum of energy of 
individual crystals greater than 20 MeV are considered for further analysis. This condition 




^WsJ^'^sl Opening angle a„Dc,sEc [deg] 

Figure 5.3: (left) Difference in the azimuthal angle ^4>mdc.psb of the exit coordinate of the 
MDC and the PSB detectors, (right) Distribution of the opening angle aMDC.SEC between cluster 
position in SEC calculated from the extrapolation of the helix from the MDC and the real cluster 
position measured in the SEC. The small spike structures seen in the spectrum are due to the 
granularity of the detector. In both figures red lines indicates the conditions under which clusters 
were assigned to one track. 
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Figure 5.4: (left) The difference in the azimuthal angle l^<t>psB,SEC of the cluster in the PSB 
and cluster in the SEC detectors. The bumpy structures on the spectrum are due to the granularity 
of the detector, (right) Difference between the time of the cluster in PSB and time of the cluster 
in SEC. In both figures red lines indicates the conditions under which clusters were assigned to 
one track. 



allows to remove low energy noise signals and as a result suppresses the background. 

The PSB clusters produced by one particle are build from neighboring detection ele- 
ments which are in a time window of 10 ns. The cluster finding routine allows maximally 
for three detection elements to build one cluster. Similarly as for the SEC the procedure 
is searching for most energetic hits, and then checks for hits in closest detection elements. 
The minimal deposited energy of a hit which is taken into account is 0.5 MeV. 

The MDC clustering and reconstruction algorithms rely on finding the hits in 17 layers 
of straw tube detectors. The procedure is realized in two steps. First the routine uses 
the pattern recognition where hits in each tube are projected into plane perpendicular to 



the beam axis 116 . Then the procedure fits a circle, and minimizes the weighted sum 
of distances between the hit position and the center of the circle. In the second stage, 
algorithm fits a straight line to hits in the 'Rz' plane creating a three dimensional helix. 
Reconstruction of particle momentum vector components is done from the curvature of the 
helix assuming the homogeneous magnetic field inside the detector. 

After obtaining cluster information from all central detectors the track finding algo- 
rithm is used. This procedure matches clusters from MDC, PSB and SEC detectors and 
assigns them to one particle. First the helix from the drift chamber is extrapolated to pass 
through PSB and SEC. Next, the procedure matches exit position of the helix from the 
MDC with the cluster position in PSB, by checking the difference in the azimuthal angle 
(j). The clusters are grouped into one track if the \lS.(j)MDC,PSB\ = \4'MDC — (I^PSb] is less 
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than 10°. The experimental spectrum of the A(f)MDC,PSB is shown in Fig. 5.3 (left). 

Furthermore, the matching between the MDC and SEC is done by checking the opening 
angle a m DC, SEC between calculated position of the hit in the calorimeter from the extra- 
polation of the helix to the SEC, and measured cluster position. In this case, the procedure 
matches the clusters into one track if the opening angle oimdc,SEC is less or equal to 20°. 
The experimental distribution of the cxmdc,SEC parameter is shown in Fig. 5.3 (right). 

Also, an additional matching criteria are checked between clusters in PSB and SEC 
detectors. In this case angular and time information for clusters are available. The absolute 
value of difference in azimuthal angle A(ppsB,SEC between clusters in PSB and SEC has to 
be less than 20°, to be accepted as originating from one particle. Also difference between 
time of the cluster registered in PSB and time of the cluster in SEC, has to be in time 
coincidence window of 80 ns (±40 ns). The experimental distribution of the azimuthal angle 
between PSB and SEC is shown in Fig. 5.4 (left), and the time difference distribution is 
shown in Fig. 5.4 (right). 



6. Identification of the pp ppr] 
reaction 



After tracks reconstruction described in the previous Section as a next step of analysis 
the identification of protons registered in the Forward Detector and reconstruction of their 
four-momentum vectors was done. Events corresponding to the pp — )• ppr] reaction were 
identified based on the missing mass technique. The main trigger in the experiment aiming 
at selection of the pp — )• ppX reaction required two hits in second layer of the Range 
Hodoscope detector, azimuthal angle agreement between FRH, FTH and FWC, and no 



hits in Veto Hodoscope. The trigger conditions were described in details in Section 5.1 



6.1 Identification of recoil protons 

In order to select protons from pp — )• ppr] reaction one has to deal with large amount 
of background reactions and random coincidences. First the reconstruction algorithm 
recognizes tracks based on signals registered in different detectors as it was described in 
Section 5.2. Then tracks which corresponds to a signals in the FPC detector are marked 
as coming from charged particle and considered for further analysis. Next, the conditions 
to reject noise signals were applied. Signals assigned to tracks are checked to satisfy the 
25 MeV minimal energy deposit condition in the whole forward detector. Further on, tracks 
are verified if they are inside geometrical acceptance of the Forward Detector which can 
detect particles only in the polar angle range of 3° — 18°. 

The particle identification is done by means of the AE — E method. In this technique 
the energy deposited by charged particle in the first layer of the Forward Range Hodoscope 
is plotted as a function of the energy deposited in whole FRH. The protons coming from 
the pp — )■ ppX reaction are seen in Fig. |6.1| (left) as two-arm band. The condition for 
proton identification is indicated on the plot as a region within the borders of the black 
lines. After these cuts, only events, with two properly identified protons were chosen for 



for further analysis (see Fig. 6.1 (right)). 



The several discontinuity in the upper proton band seen in Fig. 6.1 (left) are due to 
passive material (8 mm plexi) which is holding elements of each layer together. When a 
particle is passing through region between layers it losses some part of its energy in plexi 
material. The detection of this loss is impossible, which is revealed as a discontinuity on the 
AE — E spectrum. One can also observe a structure at low energy deposition region, where 
the most of minimum ionizing pions are giving signals. The long tail which is localized 
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Figure 6.1: (left) The particle identification plot for Forward Detector. The correlation of 
energy deposited by charged particles in the first layer of the FRH detector as a function of energy 
deposited in whole Range Hodoscope. The superimposed black line indicates cut region from which 
event were taken for further analysis, (right) Charge particle multiplicity in Forward Detector 
after particle identification. The green shaded area indicate that for the further analysis only 
events with exactly two protons reconstructed were chosen. 
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Figure 6.2: The time difference between signals from protons measured in FTH. In order to 
minimize the random coincidence, the proton candidate pair has to have time difference within 
window of ±7 ns. 
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behind the lower band originates from nuclear interactions of protons with the material 
of the detector. Also the continuous background is due to the nuclear interactions. The 
fraction of events where two protons were identified in the Forward Detector after applying 



above conditions is equal to 26% (Fig. 6.1 (right)). 

Further, in order to properly select two protons the time difference between two sub- 
sequent signals in the Forward Detector was checked. The pair with the time difference 
Atp^P2 within time window of 14 ns (±7 ns) are selected for further analysis. The corre- 
sponding distribution of Atp^pj is shown in Fig. 6.2 One can see a very steep peak around 
ns which is coming from protons and almost flat background originating from random 
coincidences. 



6.2 Identification of the 77 meson 

After reconstruction and identification of two protons in the Forward Detector orig- 
inating from pp —7- ppX reaction one can plot a missing mass distribution according to 



formula 3.3 in order to identify the production of the rj meson. The resulting spectrum of 
the missing mass is shown in Fig. |6.3| A peak originating from the production of the r] 
meson is clearly seen on the top of a continues and broad background. The shape of the 
distribution outside the r/ peak region originates form the direct pions production. The 



peak seen in Fig. 6.3 contains not only the events corresponding to the final state of the 
decay t] — )• vr+Tr^vr^ but also other charged decays of the rj meson which fulfil the condition 
of the trigger and preselection in the Central Detector. Mainly it may be the rj — >■ ir^n^j 
decay which has a branching ratio of 4.60 %. 

In order to describe the shape of the background and estimate the signal to background 
ratio on this stage of the analysis several reaction channels were simulated. The main back- 
ground constitute the pp — )• ppir^iT^ and pp — )• pp-K^ir^-K^ processes. For both reactions 
the cross section for beam energy of Eb = lA GeV are unknown. However, from the shape 
of the excitation function for these reactions one can roughly predict that the ratio of the 
cross sections is equal to about — pp^pp^+^ — ~ JOO. Moreover the simulations of the signal 
decays tj — )• tt'^tt~tt^ and background processes like: rj — )• 7r^7r~7 and r] — )• e^e~7 were 
performed. Simulated data samples were analyzed in the same way as it was done for the 
experimental data. The results of the simulations in comparison to the experimental data 



are presented in Fig. 6.3 



In order to estimate signal to background ratio simulated missing mass spectra of 
background reactions were fitted to the data, excluding the signal region from 0.52 to 
0.57 GeV, according to the formula: 

B{mm, a, 13) = a ■ fpp^ppTT+n- {fnm) + (3 ■ fpp^ppn+TT-iTO (mm), (6.1) 

where the a and /3 denotes the free parameters varied during the fit, and functions 
/pp_^pp^+j^- (mm) and fpp^ppTr+-n—TTo{mm) indicates the missing mass spectra of simulated 
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Figure 6.3: Missing mass distribution for the pp — > ppX reaction measured at beam momentum 
of pb = 2.14 GeV/c after the preselection described in this Section and identification of two protons 
in Forward Detector (black rectangles). The solid lines indicate results of the simulations of signal 
and background channels as given in the figure. The shaded area indicates the fitted sum off all 
simulated reaction channels. 

background reactions. One can see that the simulations are in a good agreement with the 
experimental data. Obtained signal to background ratio on this stage of the analysis is 
equal to 0.6. 

The method to further improve the signal to background ratio and select investigated 
decay channel will be presented in Section [7| 



7. Extraction of signal for the 

77 7r^7r^7r~ channel 

The T] meson is a particle with an average hfe time of 10~^^ s, thus it decays almost 
immediately after being produced into lighter mesons. In the WAS A facility the decay 
products are detected and identified in the Central Detector (the detail description how 



the track reconstruction and clustering algorithm work was given in section 5.3). 

Initially the tracks originating from the decay of the r] meson registered in Central 
Detector have to be correlated in time with two protons identified in the Forward Detector. 
This is done by checking the time coincidences between protons and particles in CD. The 
time of two protons in Forward Detector is determined based on signals from the FTH 
detector which is positioned 1.5 m downstream from the interaction point. Thus this time 
has to be corrected for the time-of-flight according to the equation: 

I . 
= ipTH - (7.1) 

where / denotes the distance from the interaction point to the hit position in the FTH 
detector, /3 indicates the proton velocity and tp stands for time when proton left the 
interaction region. Finally the time of particles registered in CD is compared to the average 
time of both identified protons. 

In the Central Detector the timing for the charged particles is provided by the PSB 
detector, and for the neutral particle it is readed from the SEC. Both types of particles 
are treated separately due to different distance to the interaction point and different time 
resolution of the PSB and SEC. The PSB is a plastic scintillating detector giving a very fast 
time signals with a very good accuracy. While, the organic crystals building calorimeter 
are "slower" and are characterized by worser time resolutions. For both charge and neutral 



particles the time distributions (see Fig. 7.1) shows a pronounced peak on an almost flat 



background. In case of charged particles the selection window is 14 ns wide (it 7 ns), while 
for the neutral particles it is 30 ns wide (it 15 ns). 

Furthermore, the tracks are checked to fit the geometry of the Central Detector. The 
charged particles are identified only if at least 2 axial and 3 stereo straws gives a signal 
and this is possible only if the particle passes at least five layers of the MDC. Therefore, 
the range of the polar angle for the acceptance of the charged particles is equal to 24" - 
159°. For the neutral particles, SEC acceptance enables to register signals in the range of 
20° - 169°. 

To select candidates for the r\ — )• ■k^'k~'k^ decay, the reconstructed tracks which fulfilled 
time and acceptance conditions are taken into account. Only events with exactly two 
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Figure 7.1: Distribution of differences between time of particles registered in CD and average 
time of protons in FD corrected for the time-of-flight: (left) for charged particles registered in 
PSB, (right) for neutral clusters detected in SEC. The superimposed red lines indicate the time 
window of the event acceptance, where for charge particles it is —7 ns < Atcnc l£ 7 ns and for 
neutral — 15 ns < AtcDN < 15 ns 
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Figure 7.2: (left) Multiplicity of charged particles registered in Central Detector. The green 
area indicates multiplicities equal to 2, which is taken for further analysis, (right) Experimental 
distribution of multiplicity for neutral particles detected in Central Detector after identification of 
two protons and selection of two oppositely charged particles. 
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charged particles having opposite charges, and at least two neutral clusters are saved for 
further analysis. The corresponding graph of multiplicity for charged particles is shown in 



Fig. 7.2 (left). The distribution shows that only in about 5 % of all events there are more 
than two charged particles identified. Therefore, for further analysis we take only these 
events were exactly two oppositely charged particles were reconstructed in the Central 
Detector. 

The multiplicity of neutral particles registered in SEC, after identification of protons 



in FD and two oppositely charged particles in CD is shown in Fig. 7.2 (right). From this 
distribution one can see that in most cases the two neutral particles were registered in the 
calorimeter. The events in which number of neutral clusters is greater than two is partially 
due to splitting and hadronic interactions of charged particle with the scintillator. At 
this stage of analysis we accept all the events independently of the multiplicity of neutral 
particles candidates. 



7.1 Identification of charged pions tt^ 

The selection of charged pions bases on momentum reconstruction from the curvatures 
of particle trajectories in magnetic field in the Mini Drift Chamber combined with the 
information about energy losses in the Plastic Scintillator Barrel and the Scintillating 
Electromagnetic Calorimeter {/\E — \p\ method). 

For both detectors PSB and SEC the correlation between energy deposited in the scin- 
tillating material by charged particle and absolute value of momentum can be graphed. 
In this method for the different types of particles the correlation of energy loss and mo- 
mentum differ and leave distinct bands which can be used for separation and identification 



purpose. Figure 7.3 shows corresponding identification spectra obtained from the Monte 
Carlo simulations for the r/ — t- 7r+7r~7r*'(7r'^ — t- e^e~^) decay. For both detectors four 
densely populated areas are clearly visible. For better visualization the momenta are mul- 
tiplied by sign of charge which enable to separate in the figure negatively and positively 
charged particles. 

In the calorimeter pions and electrons deposit their total kinetic energy, which at given 
momentum is larger for light electrons than for heavier pions. This allows to impose the 
identification conditions to distinguish between both type of particles as can be seen in 



Fig. 7.4 (left). The red line imposed on the plot indicates the pion identification condition 
given by the inequality: AEsec < 1-05 ■ \p\ • q — 0.06. 

To apply the same method in PSB an energy loss is normalized to the path length of 
the particle when passing through the detector. One can also see that in case of e~ and e"*" 
particles the energy loss in the whole range of the momentum is almost constant. In this 
case for a given momentum electrons deposit less energy than pions because energy loss 
decreases with increase of velocity. Experimental spectra of pions and leptons similarly as 



in SEC, reveal separate bands as it is visible in Fig. |7.3| (right) 
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Figure 7.3: Simulated distribution of the energy deposited in the calorimeter (left), and in 
plastic scintillator (right) as a function of the particles momentum and charge. The structures on 
the left part of both spectra corresponds to particles with -1 and on the right side with +1 charge 
state. The bands corresponding to pions and electrons are marked respectively. 




Figure 7.4: Experimental AE — \p\ spectra used for particle identification in Central Detector. 
Energy deposit as a function on momentum and charge: in the SEC (left), and in PSB before 
applying identification condition in SEC (right). The superimposed lines in left panel indicate 
the cut regions used to select signal events. 



In case of this analysis first selection of pions with the conditions on the correlation of 



energy and momentum in Scintillating Electromagnetic Calorimeter are applied (Fig. 7.4 
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(left)), and further events which met this requirement are checked in the Plastic Scintillator 
Barrel Fig. 7.5 This enables to reject particle misidentified in the SEC. 



^ 0.01 
O 0.009 




1.5-0.4-0.3-0.2-0.1 0.1 0.2 0.3 0.4 0.5 
Momentum x Charge [GeV/c] 



Figure 7.5: Experimental AE — \p\ spectra used for particle identification in PSB after applying 
identification condition in SEC, with superimposed lines indicating the cut regions used to select 
signal events. 



One can also see that in the experimental spectra for the PSB detector (Fig. 7.4 (right)) 
the bands are much more densely populated than vr^. This picture can be confusing 
because the cross section for the pion production is few orders of magnitude higher than 
for electrons. But one has to stress that the identification plots are made after demanding 
at least two neutral clusters in the electromagnetic calorimeter. Thus, the number of 
pions from the direct production via channel pp — )■ ppvr'^vr" drastically decreases and they 
can only pass the requirements when two fake clusters (for e.g. due to hadronic split- 
off) are wrongly identified as photon. Furthermore, requirement of at least two neutral 
particles, enhances in the selected sample fractional number of events corresponding to the 
conversion process of photons on a beryllium beam pipe. One can suppress the background 
originating from the external photon conversion based on the correlation of the distance 
between the center of the interaction region and the point of closest approach RcA of 
two helices, and the invariant mass IMBp{e~^e~) calculated assuming that the e^e~ pair 
was created in the beam pipe. The dilepton pair originating from the conversion process 
creates small values of the invariant mass and the minimal RcA value should be around 
30 mm (for particles flying perpendicular to the beam pipe) which is the radius of the 
beam pipe. The corresponding correlation between the radius of closest approach and the 



invariant mass calculated on the beam pipe is illustrated in Fig. 7.6 (left). The conversion 
and non-conversion events are very clearly separated. One can see that conversion process 
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Figure 7.6: (left) Experimental distribution of the distance Rca between interaction point and 
the reconstructed e^e~ vertex versus the invariant mass of e+e^ pair calculated assuming that 
lepton pair originates from the beam pipe, (right) Distribution of the AE — \p\ for PSB detector 
after demanding only events with RcA < 10 mm and invariant masses of e+e^ grater than 0.07 
GeV/c2. 



enhances density population in the region of low invariant masses and Rca > 30 mm. The 
electrons and pions originating from production reactions populate larger invariant masses 
and distance RcA close to 0. Furthermore, in order to check how the conversion events 
influence the identification distributions we plotted the AE— \p\ for the PSB detector under 
the condition that Rc^ is smaller than 10 mm and invariant masses e~^e~ are greater than 



0.07 GeV/c^. The resulting spectrum can be seen in Fig. 7.6 (right). After rejecting the 
conversion event region one can see that the band is almost invisible. However, it is 
important to stress that, this restriction is not used in the further analysis of the decay 



r] —7- vr^vr vr*^. 



7.2 Identification of 7 quanta and tt^ mesons 

The reconstruction of the vr^ meson relies on registration and identification of two 



gamma quanta in the Scintillating Electromagnetic Calorimeter. As it was shown in Fig. 7.2 
(right) for fraction of about 15% of events, due to splitting, more than two clusters were 
reconstructed. For these events to select pair originating from the neutral pion decay we 
apply the chi-square test according to the equation: 

^ 2 _ (JM(7,7,-) - m^o)^ 

Xij — 2 ' {' ■'^J 
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where the m^o = 134.98 MeV denotes mass of the neutral pion, IM{'yi^j) is the invariant 
mass of two photon candidates, and fi^o is the experimental mass resolution. The invariant 
mass is calculated according to the formula: 

IMi^aj) = ./iF^+r^ = 2^E^^siJ-^, (7.3) 

where and 9^.^. indicate the energy of photon candidates and opening angle between 
them, respectively. 




Figure 7.7: The invariant mass distribution for two photons registered in the calorimeter: (left) 
for the experimental data, (right) for the simulated sample. In the case of three photon candidates 



only these pairs were taken into account for which the defined in eq. 7.2 was minimal. In 
both cases a clear peak is visible at the mass of the neutral pion m^o = 134.98 MeV. The 
superimposed solid line indicates the Gaussian fit in order to determine the resolution and position 
of the invariant mass peak. The obtained resolution for both experimental and simulated data 
equals to cr w 17 MeV. 

For further analysis we accept only these two i,j clusters for which the xfj is minimal. 
Furthermore to clear the data sample and reject split-off effect of photons we have taken 
into account only pairs of clusters with opening angles greater than 40". The photon pair 
with lower values of opening angles are mainly coming from split-off processes and are 
contributing to the background. The spectrum of the invariant mass of accepted pairs 



of photons is shown in Fig. 7.7 where distribution on the left panel was obtained from 
experimental data sample and in the right panel from the simulated data. 

Experimentally obtained distribution of invariant mass of two photons shows a clear 
peak at the mass of the vr'^ meson. To compare and tune the simulation of the detector 
response we have plotted the same spectra for the simulated sample, and fitted the peak 
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region with the gauss function. Parameters obtained for both spectra are in agreement. The 
invariant mass resolution for the experiment and simulation is almost the same and equals 
to cr ~ 17MeV. Thus the Monte Carlo simulations are well tuned to the experimental 
conditions. 



7.3 Background suppression 

After selecting complete reaction chain which reads pp — t- pprj — t- pp7r~^'ir~7r^{'y^) and 
identifying all the particles in the final state using methods described in previous sections, 
one can plot distribution of the missing mass of the pp — t- ppX reaction as a function 
of the invariant mass of the decay products 7r"'"7r~7r''(77). The relevant distribution is 



shown in Fig. 7.8 A clear enhancement in the number of entries is visible at the mass of 



the r] meson. However, a large amount of background is still present in the data sample. 
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Figure 7.8: (left) The experimental distribution of the invariant mass of the 7r+7r^7r'^(77) system 
as a function of the missing mass for the pp — > ppX reaction, after particles selection and identi- 
fication introduced in previous sections, (middle) Experimental distribution of the missing mas 
for the pp — > ppX. (right) Experimental distribution of the invariant mass of the Tr+vr^Tr'^system. 



The remaining background is caused mostly by the two processes: (1) direct production 
of two pions via pp — ?■ ppvT+vr" reaction, and (2) by direct production of three pions 
via pp —7- pp7r"'"7r~7r''(77) reaction chain. For the two pion production, the previously 
described event selection and identification method turned out to be insufficient due to 
splitting of signal in the calorimeter. Nevertheless, this background contribution can be 
further suppressed using methods which will be described in this section. 

The contribution of the direct two pion production can be identified using the mo- 
mentum and energy conservation. One can look at the spectrum of the square of the 
pp — )• ppTr'^Tr~X reaction, where in case of the signal reaction X denotes missing neutral 
pion. Reconstructed missing mass of desired reaction should be around the squared mass 
of the neutral pion m q = 0.018 GeV'^/c^, and the mass for reactions where the vr" was 



not produced should be around zero. Left panel of Fig. 7.9 shows corresponding missing 
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mass spectrum for simulated signal and background reactions (solid lines), as well as the 
experimental data (black points). One can observe two peaks. First located around zero 
corresponds to the background events where no additional particle was produced (or it 
was massless). The second peak corresponds to three pion production in the final state 
located around squared vr*^ mass. To further suppress the background from the two pion 




Figure 7.9: Distribution of the missing mass squared to the system pp — > ppn^Tr^X (left) 
and missing mass to the pp — )• n^n^X reaction as obtained after the requirement that missing 
mass on the left figure is larger than 0.005 GeV ^/c^ (right). In both panels solid lines indicate 
reconstructed events from simulations of signal reaction pp — > pprj — > pptt^tt^tt'^ — > ppn^Tr^jj 
(red histogram), background 77 decay via channel pp — >■ ppij — ^ ppn'^Tr^j (blue histogram), three 
pion direct production pp — ppn'^ Tr~ n'^ — pp7r"'"7r~77 (green histogram) and two pion production 
pp — )■ ppn^T:^ (dark green histogram), black points denote the experimentally measured distribu- 
tions, and the yellow area is sum of all simulated reactions. The superimposed line in both plots 
indicates the cut chosen to suppress two pion background. The spectrum for X — ppn^ is spread 
below the sum of masses 2 • nip + w^o because of the finite experimental resolution of determining 
the four momentum vectors. 

final state, all events with squared missing mass of the ppvr+vr^ system lower then 0.005 
GeV'^/c'^ are rejected. The applied condition is indicated as a solid red line in Fig. 7.9 



(left). It reduces the 90% of the r] — t- tt^tt^j background and 97% of the pp — ?■ ppn^ir^ 
decreasing the efficiency for the signal by 41%. 

The remaining background can be suppressed by looking at the missing mass for the 
pp — )• 'ir~^Tr~X reaction, where for the signal {X = ppir^), inx should be larger than 
2 • nip + m^o = 2.01 GeV/c^. Whereas, in case of the direct production of the two charged 
pions the missing mass is only due to two protons and the distribution is shifted towards 



the lower masses. The situation is illustrated in Fig. 7.9 (right) where the experimental 
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and simulated data are shown in the same way as it was for the left panel. One can see 
that in the measured data the signal peaks are separated from each other. The cut was 
chosen to accept missing masses greater than 1.95 GeV/(P. This allows to reduce the 
remaining two-pion background to negligible level, decreasing the efficiency of the signal 
reconstruction by 0.1% only. 

The result of all previously applied condition is illustrated in Fig. |7.10| It can be seen 
that the multi-meson background (mostly direct two pion production) was significantly 



reduced compared to the one shown in Fig. 7.8 The signal to background ratio was 



improved from 0.8 to 2.2 The spectrum of the missing mass (see Fig. 7.10 (middle)) shows 
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Figure 7.10: (left) Experimental distribution of the invariant mass for the tt+tt^tt" system as 
a function of the missing mass for the pp — > ppX reaction after applying background suppression 
described in the text, (middle) Experimental missing mass for the pp ppX reaction, (right) 
Invariant mass of the Tr+Tr^Tr" system. 



a peak located at the mass of the r] meson with a tail towards the lower masses. The tail 
and continuous background under the peak is originating mostly from direct three pion 



production. The invariant mass distribution of the 7r"'"7r~7r'^ is shown in Fig. 7.10 (right). 
It also shows a maximum located around the mass of the r] meson. But the invariant mass 
distribution is much broader than the missing mass spectrum. This is because the four 
momentum vector reconstruction of particles registered in Central Detector is less accurate 
than in the Forward Detector. 



7.4 Kinematic fit 

The identification of the reaction pp — )• ppr] and the subsequent decay ry — )• vr'^Tr^vr" is 
based on the measurement of momentum vectors of all final state particles and reconstruc- 
tion of the missing and invariant masses. A finite resolution of the four momentum vector 
determination is refiected in the population of kinematically forbidden regions of the phase 
space. This obstacle will be corrected by performing a kinematic fitting of the data. In this 
experiment we have measured more variables than needed for complete description of the 
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kinematics of the final state. This redundancy will be used to improve the mass resolution 
and in consequence to improve the signal to background ratio, and to correct each event 
to follow kinematically allowed phase space region. 

The kinematic fit procedure is a least squares method which is based on variation of 
kinematic observables of all final state particles in the range of experimental resolution in 



such a way that after the change they fulfil the kinematical constraints 117 . Ax value 
is used as a measure of differences between values of a measured and the varied observ- 
ables which matches the assumed hypothesis. As a best solution fulfilling the required 
criterion one take this minimizing the function. In our case we require that event fulfil 
kinematics of the pp — )• ppir^-K^ir^ — )• pp-K^TT~^^ reaction chain. The basic conditions for 
the minimalization includes the check of the momentum and energy conservation. Also, 
conditions such as masses of intermediate particles can be introduced additionally into the 
fitting procedure. 

The x^ quantity which describes the agreement between fitted and measured variables 
reads: 



X 



2 - {Ff - Fm)S-\Ff - FMf + XcC + XuU, (7.4) 



where the Fp and Fm represent vectors of fitted and measured variables, respectively, S 
denotes the covariance matrix of the uncertainties, the Ac and stand for Lagrange 
multipliers to introduce the constraints for energy and momentum conservation C, and 
additional conditions U. 

In case of this thesis after identification of particles masses their momentum vectors 
are described by a set of three variables: kinetic energy {Ek), polar angle {9), and the 
azimuthal angle (0). For investigated reaction chain one has six particles in the final 
state: pi, p2, vr"'", tt~ , 71, 72. Therefore, taking into account the four equations for the 
momentum and energy conservation, the minimal number of independent variables which 
has to be measured to describe the system is equal to: 

A?-™" = 3 • 6 - 4 = 14. (7.5) 

However, we have measured all the variables, so we have A^^ = 18, and therefore for each 
event we have redundant information allowing us to perform kinematic fitting. The number 
of degrees of freedom can be calculated as: 

Nndf = N,- A^r" + Nu, (7.6) 

where the Njj denotes the number of additionally introduced conditions. In case of 
this analysis we introduced the condition that two photons registered in electromagnetic 
calorimeter are originating from the decay of the neutral pion. Thus according to equa- 
tion 7.6 the number of degrees of freedom reads N^df = 5. Further, one can introduce 
the probability defining the goodness of the fit for a given x^ value and a given number of 
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degrees of freedom N^df, as: 

Pf{x\ Nndf) = 1 - Fix", Nndf) = I ^ / e-H^-'Ut, (7.7) 



where F{x^ ,Nndf) is the Cumulative Distribution Function 118 , and the F denotes the 
Euler special function. The distribution of the probability of the fit Pp for the proper 
hypothesis in case of Gaussian errors should be flat 

As mentioned above to use the kinematic fit procedure one has to introduce an error for 
each measured observable. These inaccuracy depends on the particle type, its energy, polar 
and azimuthal angle, and the specific detection properties of different detectors. In order 
to estimate these errors using the WASA Monte Carlo package a sample of single particle 
events were simulated. The single particle events were chosen due to the technical reasons, 
because they ensure the uniform distribution in the whole acceptance of the detector in 
comparison to a specific reaction like for e.g. the pp — )• pprj, where particles more often are 
emitted in forward direction. The reconstruction methods and condition for the detection 
of particles are identical as in the previously described analysis. 

The convention for error parameterization for each particle was chosen as fallows: 

ae = e-"^^ - 93^'', (7.8) 

where superscripts rec and gen denote the reconstructed and generated value of the vari- 
able, respectively. The kinetic energy is parameterized by a relative difference of the 
reconstructed and generated values, and for the angles the absolute values are used. The 
errors are taken as a one standard deviation calculated from fitting the resulting distribu- 



tions given by equations 7.8 with a Gauss function. Finally all errors are determined as 
double differential functions which depends on the reconstructed kinetic energy and polar 
angle. 

The bin width for the parameterization was chosen separately for different type of 
particles. For protons in Forward Detector size of kinetic energy and polar angle intervals 
are equal to 50 MeV and 1°, respectively. For the charged pions in Central Detector the 
bin width for energy is 50 MeV and for polar scattering angle 4°. The different approach is 
for parameterization of errors of energy and angle for photons registered in the calorimeter. 
Here the energy is divided into intervals with size 50 MeV, but instead of the polar angle 
the 24 bins corresponding to number of ring building the calorimeter were taken. The 



summary of the bin size chosen for the parameterization is given in Tab. 7.1 



^It is because the probability density distribution of the cumulative distribution is always uniform 
independent of the probability density function. 
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Table 7.1: Summary of the bin size used in parameterization of errors in kinematic fit. 
The errors as a function of kinetic energy and polar angle extracted for each type of 



particle can be seen in Fig. 7.11 The plots are aligned in rows from the top for: protons, 
photons, vr"'" and tt". The columns from the left are assigned accordingly to: kinetic energy, 
polar angle and azimuthal angle. 

For protons it can be seen that relative error of kinetic energy does not depend on 
the scattering angle 6 and for energies up to 360 MeV it is almost constant (around 2%), 
for higher energies it starts to increase linearly up to 11% at 700 MeV. The increase is 
associated with the fact that protons with energy higher than 360 MeV punch through 
the FRH. The error of the polar angle is nearly independent of the energy but it increases 
with the scattering angle. Whereas, the error of the azimuthal angle decreases when polar 
angle increases, and it is almost independent of the energy. 

For photons the relative energy error decreases with increasing energy with visible 
enhancement in the regions where the central part of the calorimeter meet the backward 
or forward part. In case of the azimuthal and polar angle the errors decrease with increasing 
energy of photons. For both variables errors are increasing in the end caps of the calorimeter 
staying rather constant in the central part. 

In the case of charged pions, errors for all kinematical variables in the range from 40" to 
140° shows dependence only on the kinetic energy. For angles below 40° and above 140° the 
errors are instantly increasing. This illustrates some features of the reconstruction method 
in the MDC. These two regions reflect the front and rear part of the drift chamber. Particles 
flying through these parts fire less straws which are contributing to the track reconstruction 



and thus the inaccuracy of the measurement increases (see Fig. 7.12). 

Having the errors parameterized one can execute the kinematic fitting procedure both 
on the measured and simulated data samples. Observables in each event will be varied such 
that the momentum and energy conservation and additional constraint that two detected 
photons originate from the decay of the neutral pion are fulfilled. The obtained and 
Ppix^iNndf) spectra for simulated and experimental data are shown in Fig. 



7.13 



The distribution of the for the simulated data presented in the upper left part 
shows an significant enhancement at large values of in comparison to 



7.13 



of Fig. 

the theoretically predicted distribution for five degrees of freedom. This reflects the 
non-Gaussian error distribution which were assumed and taken into the fitting routine. 
The problem were seen for charged pions flying to the forward and backward part of the 



MDC (see Fig. 7.12). The right panel of Fig. 7.13 shows the distribution of probability 
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Figure 7.11: Distributions of the double differential parameterization of errors in determining 
the Ekin, S and (f> (from the top in rows) for: protons, photons. 7r+ and tt^. The columns from 
the left denote errors for kinetic energy, and polar and azimuthal angle, respectively. 
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Figure 7.12: Two exemplary distributions of the relative difference of the reconstructed and 
generated values of the kinetic energy of 7r+ for: (left) front region of MDC for polar angles 
24° < 9 < 28° where less straws fires and errors are no-Gaussian, (right) and middle region of 
MDC for 120° <6< 124° where errors are Gaussian. 



Pf{x^ 1 Nndf) which has an enhancement at zero corresponding to events seen in tail of the 
\^ distribution. Events which are populating low probabilities are wrongly reconstructed. 
The probability distribution is relatively flat for values of probability larger than 0.1, and 
only these events will be taken in further analysis. 



The lower two plots in Fig. 7.13 illustrate the and probability distributions obtained 
for the experimental data sample. One can see that here the Pp{x^ , Nndf) distribution is 
peaked at lower probabilities even more strongly than result obtained from the simulations. 
This is because in the data sample still a background events are present for which the fit 
hypothesis is not justified. Also the has much more enhanced tail in comparison to 
Monte Carlo spectrum due to non-Gaussian error estimate for the charged pions which 
are coming also from the background channels. To suppress majority of the background 
and wrongly reconstructed events the region of Pf{x^-, Nndf) > 0.1 was selected for further 
analysis. 

The missing mass spectrum before and after applying the kinematic fit is shown in 
Fig, 7.14 (left). One can see a slight improvement in the resolution. The ratio of the signal 
to background increased from 2.2 to 2.9. The invariant mass distribution of the tt'^-k^tt^ 
system is shown in Fig. 7.14 (right). A very narrow peak is seen in the mass of the r] 
meson in comparison to the distribution before the kinematical fit. The resolution of the 
invariant mass after the fit is equal to o" = 4 MeV/(?, compared o" = 25 MeV/c? before the 
fit. Thus, kinematical fit had improved significantly the resolution of invariant mass, and 
as a consequence also the signal to background ratio. The remaining background is due to 
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the direct three pion production. One can also see that the shape of the distribution after 
the fit is very similar to the missing mass of the two protons. This is because in the fit 
the information of protons four-momentum vectors is used, and they are determined with 
better precision than four-momentum vectors of other particle. 

Furthermore, it is also worth to notice that one can improve the resolution of fit pa- 





Probability P(z^N „) 

Figure 7.13: The (left) and probability Pp{x'^,Nndf) (right) distribution for the simulated 
(upper row) and experimental data (lower row) for the reaction pp — > pprj — > pp7r+7r^7r°(77). 
The superimposed red curve indicates the theoretical distribution calculated for the five degrees 
of freedom. The disagreement between theoretical and experimental distribution reflects the 
non-Gaussian errors taken into the fitting procedure and treated as Gaussian. 
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Figure 7.14: Missing mass distribution (left) for the pp — !■ ppX reaction and invariant mass 
spectrum for the 7r+7r^7r° system (right) before (red histogram) and after (green histogram) 
kinematic fit and rejection of events with Pf{x^, ^ndf) less than 0.1. 

rameters by demanding the invariant mass of the decay particles to be equal to the mass of 
the r] meson. However, this would result in, squeezing the whole distribution to the value 
of the mass of the r] meson, and thus subtraction of the background would be no longer 
possible. Therefore in this analysis we did not used this condition. 



7.5 Estimation of the analysis efficiency 

To determine the overall reconstruction efficiency of a particular reaction channel one 
has to take into account the geometrical acceptance of the detector and efficiency of the 
applied analysis chain. The geometrical acceptance can be studied using kinematic event 
generator by superimposing on the scattering angle of final state particles the condition 
that it has to be within a range of the detector. The geometrical acceptance was studied 
in Sec. lOl 

In order to study the efficiency of applied analysis, and the resulting background sup- 
pression, various reactions were simulated using the WASA Monte Carlo package and 
analyzed by the same software as it was used for the experimental data. The list of re- 



actions together with the number of simulated events is given in the Tab. 7.2 Second 
reaction listed, may be misidentified as a signal due to the splitting of clusters in the elec- 
tromagnetic calorimeter. Third reaction can imitate signal due to misidentification of e"*" 
and e~ as pions and splitting of signals in the calorimeter. Fourth reaction constitute a 
physical background of investigated decay. The fifth process can obscure the signal channel 
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No. 


Reaction 


Number of generated events 


1 


PP - 


ppT] 


— y pp7r~^7r TT^ — y pp7^~^7r (signal) 


20 X 10^ 


2 






PP —7- pprj —7- ppn^TT^^ 


5 X 10^ 


3 






PP —7- pprj —7- ppe^e~j 


20 X 10^ 


4 




PP 


—7- ppTT^TT^TT^ —7- pp7r+7r^77 


20 X 10^ 


5 






— )• pp-rr^TT^ 


58 X 10^ 


6 






PP — )• ppir^TT^ — )• pp^^^j 


18 X 10^ 



Table 7.2: List of simulated reaction cliannels for efficiency studies. The number in right column 
indicates the number of initially generated events. 




"J 



Figure 7.15: Reconstruction efficiency obtained after subsequent application of conditions indi- 
cated at the bottom of the figure. Results are shown for simulations for reactions listed inside the 
figure. 

by the bremsstrahlung of charged particles in the calorimeter. Sixth reaction can simulate 
a signal due to the external conversion of 7 quanta and the subsequent misidentification 
on the electron and positron as pion pair. 

The efficiency at each analysis stage can be studied by determining the number of events 
left after applying the particular condition. Figure |7.15| presents obtained reconstruction 
efficiency for several simulated reactions in the subsequent steps of the analysis. The 



name of bins corresponds to the name of applied cut and they are explained in Tab. 7.3 



Furthermore in the right column of Tab. 7.3 the efficiency of individual conditions for the 
signal reaction is given. 

As expected first biggest drop in the efficiency comes from the preselection of candi- 
dates for the investigated reaction channel in Forward and Central Detectors. As it was 
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Cut 


Cut description 


Section 


Efficiency 


JT iV-LJOl-JlJ 


r i CoCicU LlUll Ui IWU Ul U LCllo 111 1 LJ ^ LWU 

oppositely charged particles in CD and 
two OT" Tnni'P plimtPTS m SIT^((; 


Sec. 


6.2 






CDN>2 


Selection of at least two clusters 

ITl (;T~) Wltn TTl 1 Tl 1 TTl 1 1 TTl PTIPVP'V OT /T) IVTpV 


Sec. 
Fig. 


7 1 

7.2 




71.6% 


CDC±=2 


Reconstruction of tracks corresponding to 
two "navtiplps witn OTiTiomtp pnai'P'P^ 


Sec. 
Fig. 


7 1 

7.2 




47.2% 


SEC 


Particle identification in the 
SpTntillj^tTnc T^JpptT'OTnaP"'nptip (laloviTTiptp'r 


Sec. 
Fig. 


7 1 

7.4 




75.5% 


PSB 


Particle identification in Plastic 

Svr'Tnf illfif ov RfiTTpl 

kJ V'llllj llld LiVJl J_J dl 1 CI. 


Sec. 
Fig. 


7 1 

7.5 




71.9% 


27 


Selection of two photons originating from 

flip c\f^CP\\T of TT^ TTlPtlO'n 
LllC U.CL^Cl y \JL 1 \ lllCoVJll . 


Sec. 
Fig. 


7 

7.7 




92.7% 


MM2(pp7r+7r") 


Cut on the missing mass squared of the 
pp — }■ ppiT^ 'K~ X reaction. 


Sec. 
Fig. 


7 

7.9 




59.3% 


MM(7r+7r-) 


Cut on the missing mass for the reaction 

pp — )■ TT^TT' X. 


Sec. 
Fig. 


7.3 
7.9 




99.9% 


KFit 


Kinematical fit and rejection of events with 
probability Ppix^Nndf) < 0.1. 


Sec. 


7.4 




41.8% 


Total 


Final reconstruction efficiency 




0.89% 



Table 7.3: Reconstruction efficiency of the individual steps of the analysis for the signal reaction 
pp — > pprj — > 7r+7r^7r'^(77) obtained from the Monte Carlo simulations. 



shown in Sec, 4.4 the geometrical acceptance of the detection apparatus causes the lost of 
65% of all events. The rest of the decrease in the efficiency on the preselection level can 
be attributed to selective conditions in the Central Detector specially the conditions for 
identification of charged pions for which we demanded signals in all central sub-detectors. 
Second major decrease comes after applying the condition on the probability Pp(x^, Nndf) 
of the kinematical fit. 

The performed analysis allowed to increase significantly the signal to background ratio. 
The background from the two pion production was suppressed to the negligible level, and 
the signal of most dangerous decay (r/ — t- 7r^7r^7) was suppressed by more than four orders 
of magnitude. The event sample resulting after applying all conditions, contains nearly in 
100% the T] — )• vr+vr^Tr'^ decay and the only background left corresponds to the same final 
state originating from direct three pion production. This background can be seen in the 
missing and invariant mass spectra shown in Fig. |7.14 as a continues distribution spreading 
over much larger range than the rj meson peak. In the next chapter we will introduce the 
method how to estimate contribution of this background using a fit of polynomial function 
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to the missing mass distribution. 



8. Results for the 77 tt^tt^tt decay 

Identification of all final state particles of the rj — )• 7r~^TT~iT^ decay enables to use a 
Dalitz plot analysis to study the dynamics of three body system which manifests itself 
in the population density distribution in this plot. The parameterization of the events is 
given conveniently in the normalized variables X and Y, dependent on the kinetic energies 



of three particles in the rest frame of the rj meson which were introduced in equations 2.8 



and 2.9 In this chapter we will present the final results for the Dalitz plot distributions. 



8.1 Resolution of the observables determination 

In order to use the Dalitz distribution one has to know the resolution for the deter- 
mination of the X and Y observable and decide about the bin width. The interval size 
depends on the collected statistics and the accuracy in determination of four momentum 
vectors of particles in the final state. 

First to check the agreement between simulations and experiment we have determined 
distribution of the difference between the values of X and Y before and after the kinematical 
fit for simulations and measured data: 



AX = XkFU — Xrec, (8.1) 
Ay = YKFtt - Yrec. (8.2) 



The result is presented in Fig. 8.1 where the left panel shows the distribution of AX and 
right of Ay. The simulated data were normalized to the experimental distributions such 
that integrals of both are the same. It is clearly seen that for both observables a good 
agreement between simulated and measured data were achieved. However, one can also 
notice that the distribution of Ay variable is not symmetric and shifted to lower values. 
One of the reasons of the asymmetry in the Ay variable could be high and rapidly changing 
values of errors at the edges of the acceptance for each type of particle as it was shown in 



Fig. 7.11, and the lack of the availability in present WASA-at-COSY software of inclusion 
into the fit of a full covariance matrix. Observed effect influence the tiny values of searched 
asymmetry parameters. Therefore, for the further extraction of these parameters in the 
r] — )• vr^vr^vr'^ decay, we will use values of X and Y as reconstructed before the fit, which 
are independent of the simulations of errors, and their correlations In order to, avoid 

^It is worth to mention that the kinematically fitted events for the rj — > tt^tv'tt'^ will be further used 
in the analysis of the t] — > yr^e+e^ where the three pion decay is used as a normalization channel for the 
estimations of the branching ratio. In this case kinematic fit improves significantly the resolution of the 
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KFit " rec ' KFit " rec 



Figure 8.1: Distribution of differences for (left) X and (right) Y variable between values after 
and before the kinematical fit. Red histogram stands for the experimental data and green for the 
Monte Carlo simulations. The simulated distributions were normalized to the experimental data. 



the large uncertainties of the energy and polar angle for charged pions, at the edges of the 
acceptance, for the determination of the asymmetry parameters, we accept only pions with 
kinetic energy greater than 50 MeV and with scattering angle greater than 35°. 

In order to estimate the resolution of the reconstruction of the X and Y we have taken 
simulated data, and for each event compared generated values with values reconstructed 
after the kinematical fitting: 

^Xres = XkFU — ^gen, (8-3) 



AYres = YkFU - Ygen, (8.4) 

where the subscript gen stands for the value calculated from generated four momentum 
vectors, and KFit denotes the values calculated from momentum vectors reconstructed 
after the kinematical fit. The resulting distributions of AX res and AYres were fitted with 



the Gauss functions. The procedure is illustrated with the spectra shown in Fig. 8.2 
Obtained resolution for the determination of X is equal to ax = 0.061, and for Y is 
(Ty = 0.055. Therefore, the binning of the Dalitz plot histogram was chosen to be in 
intervals of 0.2, which are roughly equal to about three standard deviations. 



invariant mass distribution, and and in turn it improves also the identification of the decay channel and 
allows for the more effective background suppression. 



8.2. Dalitz Plot 
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■KFit " ■^true ' KFit " 'true 



Figure 8.2: The distributions of AX^es (left) and AYres (right) obtained from simulations. 
The superimposed Hues indicate the result of fit with the Gaussian function. 

8.2 Dalitz Plot 

The final data sample used to determine the Dalitz plot contains predominantly only 
events from searched pp — t- ppr] — t- pp-K^ir'^'ir' process. The background which is left in 
majority originates from the direct production of three pions. In order to plot the Dalitz 
distribution this background has to be subtracted bin by bin. 

To remove the background one could use the Monte Carlo simulations, but for beam 
energy in range of 1.4 GeV neither total nor differential cross sections for the three pion 
production are known, thus it is not possible to describe the mass distribution of the 
vr^vr^vr'^ in a reliable way. Therefore, we will subtract the background using the method of 
polynomial fit outside the peak region. We apply the formula of a fourth order polynomial 
function, which reads [119]: 



f{mm, a,b,c,d,e) = a + b ■ mm + c ■ mm^ + d ■ mm^ + e • mm^ (8-5) 



where a, 6, c, d, e are free parameters varied during the fit. 
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Figure 8.3: Experimental distributions of the missing mass determined for three exemplary bins 
of the Dalitz plot: (upper panel) before the background subtraction, and (lower panel) after 
the background subtraction using a polynomial fit. The fitted function is shown as a red solid line. 
Vertical dashed lines indicates the region in which the number of signal events were counted. 



The upper panel of Fig 8.3 presents exemplary missing mass for three bins of the Dalitz 
plot. The superimposed red lines denote the result of the fit with the function defined in 
Eq. |8.5| The polynomial was fitted to the data outside the maximum corresponding to 



the signal. For each presented in Fig. 8.3 missing mass spectrum, a very clear peak can 



be seen at the mass of the r] meson, above a continues almost flat background. One can 
also see that the fitted function is in a good agreement with the shape of the background. 
After subtraction of background we counted the number of events in the intervals marked 
by the dashed lines. 

Figure 8.4 (left) shows measured Dalitz plot for the rj — t- vr^vr^vr^ decay before back- 
ground subtraction. One can see that part of the events is outside the kinematical bound- 
aries of the Dalitz plot. Therefore for further analysis we will consider only event in bins 
which centers are fully inside the kinematical boundaries. The distribution on the right 
shows the result obtained after the background subtraction, performed separately for each 
bin. 

In the next step the experimental distributions have to be corrected for the detection 
and reconstruction efficiency to compare to the model expectations. The overall efficiency 



was determined in Sec 7.5 and it was defined as a fraction of number initially simulated 



8.2. Dalitz Plot 
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Figure 8.4: Experimental distributions of the Dalitz plot: before background subtraction (left), 
and after background subtraction (right). The black solid line indicates the kinematical boundary 
of the Dalitz plot. In the right panel number of events in each interval is indicated. 



events to the number of reconstructed events after the apphed analysis. In order to be able 
to draw a conclusion as much model independent as possible the efficiency correction has to 
be determined separately for each interval of the distribution of interest. Therefore, for the 
Dalitz plot the correction has been done using the two-dimensional efficiency distribution 
derived for each bin. This method enables to be independent of the model used for the 
description of the decay amplitude in the simulations. The only dependence may occur 
due to the finite bin size, but anyhow this influence is only due to the non-linearity of 
the distribution within the bin limits and this in our case is negligible. To determine the 
efficiency we have simulated a sample of events assigning a weight to each event according to 



the ChPT Leading Order predictions 34,35] (see Tab. 2.1). The relation for the efficiency 
is given by: 

J2 ■ w''-'"'{X,Y) 

where w^^'^ and iff^" denotes the weights of the events and the summation is done over 
reconstructed and generated events within given X,Y interval of the Dalitz plot. 



The Dalitz plot for model generated events is presented in Fig. |8.5| (left). Generated 
four momentum vectors of all particles were used in the Monte Carlo simulations of the 
response of the detector. Next, all events were subject to the same analysis chain as the 



experimental data. The resulting spectrum for reconstructed events is shown in Fig. 8.5 



(middle). The two dimensional efficiency determined as a function of each bin of the Dalitz 
plot is shown in the right panel of Fig. |8.5| The maximum efficiency was found for values 



of Y close to -0.5 and X values close to 0. Figure 8.6 shows the efficiency corrected Dalitz 
plot for the r] — t- tt^tt^tt^ system. The population is decreasing with increasing values of 
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Figure 8.5: Dalitz plot distribution for: (left) all generated events weighted with the ChPT 
Leading Order predictions, (middle) reconstructed events after passing whole analysis chain, 
(right) The efficiency determined as a function of the Dalitz plot. 




Figure 8.6: Final Dalitz plot after background subtraction and efficiency correction. 



Y variable which is proportional to the energy of neutral pion. 

The Dalitz plot density population is expected to be described by transition amplitude 
|Mp given by Eq. 2.10 Having the background free and acceptance corrected Dalitz plot 



dof 


a 


b 


c 


55 


-1.039 (fixed) 


0.27 (fixed) 


0.05 ±0.10 



Table 8.1: Results of the two dimensional fit to the Dalitz plot distribution of the amplitude 
|M|2 = ^2(1 + a-Y + b-Y^ + c- X). 



8.3. Determination of asymmetries 
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one can fit the two dimensional distribution by minimizing the x function: 

where iV stands for the number of the events corresponding to r/ — )• ■k^'k~'k^ decay after 
the background subtraction, and is the uncertainty of the acceptance corrected number 
of signal events. The bins indicated in white are not included in the fit. We fitted |Mp 
including terms up to squared X, however we have fixed a and h parameters as predicted by 



ChPT Leading Order. Obtained results for fit is given in Tab. 8.1 The coefficient c which 
reflects the degree of C-invariance violation is consistent with zero. As it will be shown in 
the next section this is consistent with the values of the obtained left-right asymmetry. 



8.3 Determination of asymmetries 

The integrated asymmetries of the Dalitz plot are sensitive for the C violation in 
the amplitudes of a given 1 state, therefore it is important to estimate their values. The 



definition of all three asymmetries were given in Chapter 2 by equations 2.11 2.12 2.13 



The left-right asymmetry Ai^ji is sensitive to the C violation averaged over all / states, 
however the quadrant asymmetry Aq and sextant asymmetry As can test C invariance for 
specific 7 = 1 and / = 0, 2 isospin states of the r\ — )• vr^vr^vr'^ decay, respectively. 

To determine the asymmetry parameters the Dalitz plot was divided into regions as it 



is shown in Fig. 2.2 Then according to formulas 2.11 2.12 2.13 we summed the events 



separately for odd and even regions of the Dalitz plot. Next for each summed region 
we reconstruct the missing mass of the — J- ppr/ reaction as it is shown in Fig. 8.7 
Furthermore, to determine the number of events corresponding to the r] — )• vr^Tr'Tr'' decay 
in each region we have subtracted the background using the polynomial fit method, the 



same which was applied in previous section according to the formula 8.5 Figure 8.7 
shows the corresponding missing mass distributions. Moreover, inside of each plot the 
corresponding spectrum of the missing mass after the subtraction of the background is 
shown. The number of r/ — )• 7r"'"7r~7r^ events were counted in the region denoted by dashed 



lines in Fig. 8.7 The reconstruction efficiency was estimated based on the simulations 



of signal reaction. The efficiency was estimated separately for each region of the Dalitz 
plot as ratio of number of reconstructed and generated events. Obtained values of the 



reconstruction efficiency for different regions of the Dalitz plot are collected in Tab. 8.2 



One can see that in all regions the efficiency is almost the same and it is equal to around 
0.67% 

The final results for the asymmetry parameters are: 



Alr = (+0.33 ± 0.38stat) X 10" 
Aq = (-0.18±0.38,t„t) X 10- 



(8.9) 
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Dalitz plot region 


Efficiency e 




0.6656% 




0.6659% 


ei3 


0.6719% 


£24 


0.6721% 


€135 


0.6731% 


£246 


0.6709% 



Table 8.2: Efficiency for tlie different regions of tlie Dalitz plot calculated from the simulations 
of the pp -> ppi] — > ppir^TT^TT^ weighted according to the ChPT Leading order predictions. 



As = (+0.06 ± 0.39stat) X 10-2 (8.10) 
where the second number denotes the statistical errors. The statistical error was calculated 



according to the formula for the error propagation which reads: 

2 



N%.a^{NL) + Nl.a^iNR), 



Ul) 



{Nr + NlY 

where the cr{Ni/ji) = ^/Ni^/r + Bf^/R, the Ni^jr denotes the number of signal events, and 
B^jR stands for the number of background events under the r] meson peak. The evaluation 
method of the systematic error will be given in the next section. The formula given by 



equation 8.11| can be also generalized for the calculation of the statistical uncertainties 



a{AQ) and (t{As). In right panel of Fig. 8.8 the result of this work is compared to previously 
determined values of asymmetries. The errors shown are the statistical and systematic for 
all presented measurements. 



8.4 Estimation of the systematic uncertainty 

Beside the statistical uncertainty of the measured observables one has to investigate the 
systematic effects which may change the result. To evaluate the systematic uncertainty of 

as 



established values of asymmetry parameters, we have used the method described in 120 
the method which is recommended by the WASA-at-COSY collaboration. In this work 120| 
the systematic uncertainty /S.sys is calculated as a difference between the result obtained 
with cuts used in the analysis and a result obtained after changes of these cuts. Next ^sys 
is compared to agys defined as: 

^changed ^ original^ (8.12) 

where the (Jchanged denotes the statistical uncertainty of asymmetry obtained after the 
change of the cut, and CForiginal indicates uncertainty of the originally obtained result. The 
ratio of is than used as a measure of the significance of the determined systematic 
error. If this ratio is smaller than one the systematic error due to the tested effect is 
negligible. 



8 






Figure 8.7: Missing mass distribution for the pp — pprj reaction determined to calculate the 
asymmetry parameters of the Dalitz plot. From the top in rows: (1) left-right asymmetry. (2) 
quadrant asymmetry, (3) sextant asymmetry. The red solid line indicate the fit of the polynomial 
function to subtract the background. 
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Figure 8.8: (left panel) Comparison of values of asymmetries obtained in this work with results 
determined by previous experiments, and a value given by PDG. (right panel) Number of events 
reconstructed which were used for calculation of asymmetries according to formulas 12.111 12.121 
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Below we will discuss the applied changes of the selection criteria. Each change was 
done separately for individual cuts assuming that the sources of the systematic errors are 
uncorrelated. After the change of one condition at a time the others were kept as in 
the original analysis, and the whole analysis was repeated. The obtained changes of the 



asymmetries (A A) due to variation of the condition listed below are presented in Tab. 8.3 



1. Identification of charged particles in SEC 

The identification of charged particles is based on the correlation between energy de- 
posit in the electromagnetic calorimeter and the momentum determined from the sig- 
nals measured by drift chamber. The corresponding distribution is shown in Fig. 7.4. 
The pair of charged particles is treated as a vr^ and vr" if both of them fall into re- 
gion which is below the cut lines. The cut lines are parameterized with the following 
function: AEsec = a ■ ■ q + b. To estimate the systematic influence caused by the 
charged pion identification in the calorimeter, the red lines shown in Fig. 7.4 separat- 
ing leptons from pions are moved closer to the pion bands by changing the parameter 
b. We decreased the allowed area for pions by 5%, and repeated the analysis. 

2. Identification of charged particles in PSB 

Second method of particles identification for pions is based on the energy deposited 
in the Plastic Scintillator Barrel as a function of the reconstructed momentum. The 
method is shown in Fig. 7.5. In order to estimate the systematic effect, the area 
marked on the figure by the red lines was decreased by 5%, and analysis was repeated. 

3. Suppression of the 2tt background via missing mass squared for the pp — 
ppir^TT^X reaction 

This cut was dedicated to reduce the large amount of the two pion background 
originating from the direct production. The cut can be seen in Fig. 7.9 (left). To 
estimate the influence of this cut on the final results we repeated the analysis reducing 
the condition by 0.0025 GeV^/c^ towards the lower missing masses. 

4. Suppression of the 2tt background via missing mass for the pp — 7r+7r~X 
reaction 

In the original analysis to suppress the rest of the two pion background we have used 
the condition which was presented in Fig. 7.9 (right). To establish it influence on the 
obtained result we have change the cut by 0.05 GeV/c^ toward the higher missing 
masses, and repeated the analysis. 

5. Background subtraction 

The direct production of pions obscures the interesting signal if one wants to look 
at missing mass distributions. In order, to subtract this background a polynomial 
function was fitted to the continues background outside the signal peak region. In 
order to check how the choice of the order of the polynomial changes the results we 
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have fitted the background with fifth-order polynomial function instead of originally 
used fourth-order. 



The results of performed systematic checks are gathered in Tab. 8.3 where the val- 
ues of the change of asymmetry parameters obtained after the change of individual 
cuts together with the estimated asys for the cut are given. According to the method- 
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a 
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AA 




sys 


Oryg. 


Alr 


0.0033 


0.0038 
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values 


Aq 


0.0018 


0.0038 


— 


— 
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As 


0.0006 


0.0039 
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c 


0.05 


0.10 


— 


— 
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1. 


Alr 


0.00392 


0.00401 


0.00062 


0.00129 


0.48 




Aq 


-0.00241 


0.00401 


0.00070 


0.00128 


0.56 




As 


0.00101 


0.00399 


0.00041 


0.00084 


0.49 




c 


0.072 


0.120 


0.023 


0.067 


0.34 


2. 


Alr 


0.00353 


0.00391 


0.00023 


0.00093 


0.25 




Aq 


-0.00151 


0.00392 


0.00029 


0.00094 


0.31 




As 


0.00081 


0.00380 


0.00021 


0.00087 


0.24 




c 


0.055 


0.110 


0.005 


0.04701 


0.11 


3. 


Alr 


0.00264 


0.00361 


0.00067 


0.00118 


0.56 




Aq 


-0.00101 


0.00362 


0.00079 


0.00117 


0.67 




As 


0.00110 


0.00379 


0.00050 


0.00092 


0.54 




c 


0.049 


0.099 


0.001 


0.006 


0.14 


4. 


Alr 


0.00370 


0.00378 


0.00040 


0.00039 


1.03 




Aq 


-0.00231 


0.00377 


0.00051 


0.00048 


1.08 




As 


0.00121 


0.00386 


0.00061 


0.00056 


1.10 




c 


0.046 


0.105 


0.004 


0.034 


0.13 


5. 


Alr 


0.00412 


0.00388 


0.00082 


0.00078 


1.04 




Aq 


0.00274 


0.00391 


0.00094 


0.00092 


1.02 




As 


0.00140 


0.00382 


0.00080 


0.00079 


1.01 




c 


0.060 


0.116 


0.00941 


0.05807 


0.16 



Table 8.3: Results of the studies of systematic effects. First column indicates the number of the 
investigated effect as listed in the text. Third and fourth columns shows the asymmetry value 
and its statistical uncertainty as determined in this work. Fifth column indicates the estimated 
contribution to the systematic effect (A. A), and sixth column shows its uncertainty. 



ology presented in reference 120| we consider contributions to the systematic error as 



non-significant and hence negligible if is less than one. Hence, the final systematic 
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uncertainties amount to: 

(TsysiALR) = 0.09 • 10-2, 

asysiAg) = 0.11 • 10-2, 

and systematic error for the determination of the parameter c is neghgible. For the com- 
parison we have also estimated the systematic errors more conservatively by adding in 
quadrature all contributions to the systematic error. In this case errors amount to: 

asys{ALR) = 0.13 ■ 10-^ 

asysiAg) = 0.15 • 10-2, 

asys{As) = 0.12-W~^, 

and are still much lower than the statistical uncertainty. One also has to note, that beside 
presented above sources of the systematic effects, there are still another checks which may 
be done in further studies of investigated reaction. 



8.5 Discussion 

In the first part of this thesis the charge conjugation has been tested in jy — )• Tr+Tr-vr'^ 
decay. As a main result we have obtained three asymmetries: A^ji, Aa and ^4^ of the 
Dalitz plot sensitive for the C symmetry breaking in different isospin states of the final 
particles. The estimated values of the asymmetries are listed below together with the 
statistical and systematic uncertainties: 

Alr = (+0.33 ± 02Sstat ± Omsys) X 10-2, 
Aq = (-0.18 ± 0.38stat ± 0.11,2,,) X 10-2, 
As = (+0.06 ± 0.39stat ± O.lO^j,,) x 10-2. 
One can see that obtained values of all studied parameters are consistent with zero and 



with the previous measurements (see Fig 8.8). Therefore, we can conclude that the charge 



conjugation is invariant in strong interaction at the level of achieved accuracy. 

Moreover, we have fitted the Dalitz plot density according to the phenomenological 



parameterization of the transition amplitude given by equation 2.10| As expected from 
C-invariance the c coefficient is equal to zero within uncertainty: 



C = 0m±0.10stat- 
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9. Extraction of rj 



decay 



In order to investigate the charge conjugation invariance in the electromagnetic inter- 
actions one can study the decay of the rj meson into the ■iT^e~^e~ system. This rare process 
till now has been not observed in any experiment. This decay has a similar final state as the 
decay rj — )• ■k'^-k~tt^ which was described in previous chapter. Thus, the reconstruction, 
preselection of tracks and proton identification for both processes are the same. These 
stages of the analysis were shown in chapters 5 and 6. Therefore, in next sections we will 
describe only the selection aiming at identification of the final state with two oppositely 
charged leptons e"*" and e~ and two photons originating from the vr*^ decay. 

The search for rare decays requires a good understanding of a physical processes, de- 
tector response, reconstruction methods and background subtraction. Therefore, for the 
purpose of the investigation of the rj — )• 7r^e'^e~ decay we have simulated using Monte Carlo 
methods several processes which may contribute as a background to the searched decay. 



Table 9.1 shows the simulated processes together with the number of generated events used 
in the simulations. The table is divided into three parts indicating (i) signal reaction, (ii) 
background originating from the decay of the r] meson, and (iii) background from the direct 
multi-pion production. Reactions two and nine possess similar final state as signal process 



No. 


Reaction 


Number of generated events 


1 


PP ~^ PPV ~^ ppe^e~7r^ — >• ppe~^e~^j (signal) 


5 X 10^ 


2 


pp —7- pprj —7- ppiT~^TT~1T^ —7- ppH^TT^^^ 


20 X 10^ 


3 


pp —7- pprj —7- ppTT^TT^TT^ — )• pp7r'^7r~e'*'e^7 


1 X 10^ 


4 


pp — )• pprj — )• ppTT~^TT~'y 


5 X 10^ 


5 


pp — )• ppT] — )• ppe~^e~j 


20 X 10*5 


6 


pp — )• ppr] — )• ppir^'yy — )• ppjjj^ 


3 X 10^ 


7 


pp — )• ppT] — )• ppir^'yy — )• ppe"'"e~777 


3 X 10^ 


8 


pp — pprj — PP77 


6 X 10^ 


9 


pp —7- ppTT^TT^TT^ —7- pp7r+7r^77 


20 X 10^ 


10 


pp — )• ppn^TT^ 


58 X 10^ 


11 


pp — )• pp-K^TT^ — )• pp^^^^ 


18 X 10^ 


12 


pp —7- ppir^TT^ — )• ppe'''e~777 


8 X 10^ 



Table 9.1: List of simulated reactions which may contribute to the background in the search for 
the 77 — > e^e^TT*^ decay. Right cokimn indicate the numbers of initially generated events. 



with two charge particles and two gamma quanta in the final state. The misidentification 
of the reaction could be due to mistaken identification of the pions as electrons. In the case 
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of third reaction the sphtting of signals in the calorimeter and misidentification of pions 
as electrons, and bending of electrons inside the beam pipe may be wrongly recognized 
as a signal. The fourth, fifth and tenth reaction can be misidentified as a signal due to 
bremsstrahlung effect in the calorimeter. Reactions six, eight and eleven have two or more 
photons in the final state, and these can cause the external conversion of gamma quanta 
on the beam pipe, resulting in production of electron-positron pair. In the reaction seven 
and twelve merging of clusters in the calorimeter may lead to wrong identification of these 
reactions as a signal. 

In the next section we will describe the selection criteria, established based on the sim- 
ulations, aming at suppression of background and the identification of the signal reaction. 



9.1 Identification of meson and leptons e 



The final state of r/ — )■ TT^e~^e~ decay consists from two photons originating from the 
decay of the neutral pion and two oppositely charged leptons. As a first step in identifica- 
tion of the desired decay one has to recognize all final state particles. The reconstruction 
of the TT^ meson is based on the identification of two photons in the calorimeter. From the 
analysis we accept only these events for which two or more clusters were reconstructed. 
However, part of identified clusters are due to splitting and bremsstrahlung of the charged 
particles. Thus, from the data sample one has to choose only events which contain two 
gamma quanta originating from the decay of the vr*^ meson. To this end, we will use the 
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Figure 9.1: Experimental invariant mass spectrum for two gamma quanta registered in the 
calorimeter. In the case of three photon candidates only these pairs were taken into account for 



which the defined in eq. 7.2 was minimal. The superimposed vertical solid lines at 90 MeV/c^ 



and 180 MeV/c^ indicate invariant mass region accepted for the further analysis. 
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chi-square test according to the equation 7.2 and calculate for each 7j7j pair an invariant 
mass using formula 7.3 For the further analysis we will accept only these photon candi- 
dates for which the xfj is minimal. To suppress splitting effects we restricted the opening 
angle between two clusters to be greater than 20°. The distribution of the invariant mass 



of accepted photon pairs is shown in Fig. 9.1 Furthermore, we have restricted the values 
of the invariant mass of two photons to be in the range from 90 MeV/c^ to 180 MeV/c^. 



This condition is indicated as a red solid line in Fig. 9.1 



Before the identification of electrons and positrons one has to deal with a large back- 
ground coming from the other charged particles, mostly pions originating from direct pro- 
duction reactions and the r] meson decays. In order to reduce this background we will 
use the invariant mass distribution for two oppositely charged particles detected in central 
detector. In case of e^e~ pair the invariant mass due to small mass of the electron and 
positron, the fact that they originate from the virtual photon, and the form factor mass 



dependence, the spectrum should be peaked near zero 121,122]. While for the pions the 
distribution should be shifted towards much higher invariant masses. The situation is il- 



lustrated in Fig. 9.2 where on the left panel simulations are shown for signal {rj 



) 



and background reactions. It is worth to notice that in Fig. 9.2 we present only channels 
with pions in the final state, as this cut are used to suppress multi-pion background. In the 



right panel of Fig. 9.2 the corresponding experimental distribution of the invariant mass 
is presented. To suppress the pion background we accept only these events for which the 




Figure 9.2: Invariant mass distributions of two oppositely charged particles registered in central 
detector, (left) Distribution of simulated events for rj — > 7r°e+e^ signal decay (solid black line) 
and background reactions with pions as indicated inside the figure, (right) The experimental 
distribution of invariant mass. The superimposed vertical solid line indicates the cut used for 
accepting only events with invariant masses smaller than 40 MeV/c^. 



invariant mass of two charge particles is smaller than 40 MeV/c^. The applied condition 
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is indicated as a solid red line in Fig. 9.2 It reduces the pion background almost by the 



i, decreasing the efficiency of the signal reconstruction by 24% only. 



After reduction of the pion background one can process with the identification of elec- 
trons and positrons emerged into the electromagnetic calorimeter using the AE—\p\ method 



(the detailed description of this method was given in Sec, 7.1). Charged particles scattered 



into the calorimeter deposit their total kinetic energy which at given momentum is larger 



for electrons than pions. This situation can be observed in the left panel of Fig, 7.3 where 
deposited energy in the calorimeter by charged particles as a function of their absolute 
momentum for simulated 7] — )■ 7r~'"7r~7r'^(7r'^ — )• e^e~7) decay is shown. The four clearly 
separated bands can be seen. The experimental distribution after applying previous cuts 



is shown in Fig. 9.3 In this figure only bands coming from electrons and positrons are 
clearly visible since pions are strongly reduced due to application of cut on the invariant 



mass of charged particles described previously (see Fig. 9.2 ). However, still a small fraction 



of pions can be seen in the regions bellow the e"*" and bands. Therefore, in order to 
improve the selection of electrons and positrons, we imposed an identification condition 
given by an inequality: AEsec > 1-05 - {pi ■ q — 0.06, which is indicated as a solid red line 
in Fig. 9.3 




Momentum x Charge [GeV/c] 



Figure 9.3: Experimental AE — \p\ spectrum used for electron and positron identification. The 
AE corresponds to energy deposited in the electromagnetic calorimeter and the \p\ was recon- 
structed based on tlie bending of the trajectory in tlie magnetic field. The superimposed lines 
indicate a cut used to select candidates for electrons and positrons. 



9.2. Suppression of split-ofF and conversion background 



9.2 Suppression of split-off and conversion 
background 
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Several reactions can obscure signal of searched decay due to split-off effect caused by 
the Bremsstrahlung of a charge particles in the calorimeter. The bremsstrahlung results in 
emission of an additional photon when a charge particle is passing close to the nuclei. In 
such a case the resulting photon emerges with a small opening angle to the parent particle 
and with small kinetic energy. Therefore, for example in case of the decay r] — )• e~^e~'j 
one additional gamma quantum arising from the bremsstrahlung can imitate the signal 
channel. In order to suppress this type of background we will restrict the smallest invariant 
mass of a charged and neutral particle pairs to the values larger than 120 MeV/c^. The 
smallest invariant mass was taken as a lowest value of the mass calculated for all the 
combinations of pairs of charged and neutral particles candidates. In case of the signal 
reaction (r/ — )• Tr^e~^e~ — )• ^je~^e~) an average relative momentum between 7 quanta and 
leptons are much larger than it is in the case of wrongly identified splitted clusters. The 
corresponding spectrum of invariant mass obtained from simulations for signal and several 
background channels is shown in Fig. |9.4| (left). Results of simulations are shown for these 
reactions which may be misidentified as signal due to the splitting effects in electromagnetic 
calorimeter. One can see that in case of the split-off reaction the invariant mass is peaked 
close to zero due to small opening angle and small energy of photon relative to the charged 
particle. For the signal decay the invariant mass distribution is much broader and shifted 
towards higher invariant masses. The experimental distribution of the minimal invariant 




Figure 9.4: The smallest invariant mass of pairs of charged and neutral particle candidates: 
(left) obtained from the simulation, and (right) from the experimental sample. The solid red line 
denotes the cut used in order to suppress events with the bremsstrahlung effect. 
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Figure 9.5: (left) Illustration of a primary vertex originating e.g from the 77 — > e+e^7 decay 
occurring in the interaction point (IP), (right) Diagram of an external conversion of a photon in 
the material of a beam pipe resulting in production of an electron-positron pair. 



mass of a charged and neutral particles is presented in Fig. 9.4 (right). To suppress events 
with bremsstrahlung photons we accept only events for which the invariant mass is larger 
than 120 MeV/c^. 

Another important effect which has to be studied when considering rare decays of the 
rj meson containing electrons is the external conversion of photons in the detector material 
(see also Sec. 7.1). The 7 quanta emerged from the interaction point before reaching the 
electromagnetic calorimeter have to pass through a Beryllium beam pipe with a thickness 
of 1.2 mm. However, when a photon passes through the material it can interact with the 
nuclei, and convert into e~^e~ pair, which in specific cases can obscure the signal from the 
r] — )• TT^e'^e~ reaction. In case of electron and positron emerging from the interaction 
point (IP) their momentum vectors calculated under the assumption that they originate 
from the beam pipe, will have "wrong" directions due to the curvature of the tracks in 
the magnetic field (see Fig. 9.5 (left)). As a result, the invariant mass calculated based 
on these vectors will be larger than the mass of e~^e~ pair emitted from the beam pipe 
region. In contrary, for the case of e~^e~ pair originating form the 7 conversion at the 
beam pipe the momentum vectors at the beam pipe will be almost parallel to each other 



(see Fig. 9.5 (right)), thus resulting invariant mass will be relatively small. Therefore, to 
suppress the background from pair creation one can use the correlation of the distance 
between center of the interaction region and the point of the closest approach of two 
helices reconstructed {Rqa) and the invariant mass of a lepton pair calculated under the 
assumption that they were created in the beam pipe. The corresponding experimental 
distribution of the RcA as a function of the invariant mass of the e~^e~ pair is presented 



in Fig. 9.6 One can see two regions populated on the plot. One which is located at small 
values of invariant mass and radius above 30 mm, and second around invariant mass values 
of 20 MeV/c^ and radius below 10 mm. First mentioned region corresponds to the external 
conversion of 7 quanta and, the second region to the reactions where rj meson decayed in 
the interaction region (IP) into channel with e~^e~ pair in the final state. However, the 



9.3. Further background reduction 



93 




Figure 9.6: Experimental distribution of the distance Rca between the interaction point and 
reconstructed e^e^ vertex as a function of invariant mass of lepton pair calculated assuming that 
this pair originates from the beam pipe. The red solid line indicate the cut applied in order to 
suppress the external conversion. 



conversion region is not so enhanced as it was seen in Sec. 7.1 This is because on this 
stage of the analysis we have already strongly suppressed events with photons which might 
cause the conversion by restriction on the minimal invariant mass on charged particle and 
neutral particle candidate. To suppress the conversion background we applied a cut which 
is indicated as a solid red line in Fig. 9.6 The condition results in rejection of nearly 100% 
of all conversion events. 



9.3 Further background reduction 

At this stage of the analysis one has complete reaction chain identified: pp — )• pprj — t- 
e~^e~iT^ —7- e'''e~77 with all the particles in the final state. However, to extract events 
corresponding to the interesting decay one has to further suppress background originating 
from several reactions. For that purpose one can use the missing mass of two protons and 
the invariant mass of decay products. In this section we will introduce further selection 
criteria aiming at the background suppression. 

First we will strongly restrict invariant mass of two photons identified as originating 
from the neutral pion. This condition will allow to reduce the remaining background 
from the misidentified reactions like: r] — )• tt^j^ — )• e'^e~777 with additional gamma 



quanta. Figure 9.7 (left) presents the invariant mass of two registered photons for simulated 



background and signal processes. To limit this background we accept only masses in the 
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Figure 9.7: Invariant mass of two gamma quanta for: (left) simulations, and (right) measured 
data. The red solid line indicate region of invariant masses from 120 MeV/c^ to 150 MeV/c^ used 
to select candidates for the signal reaction. 



region from 120 MeV/c^ to 150 MeV/c^ which corresponds to the FWHM of the simulated 
signal distribution. The spectrum of the invariant mass from experimental data can be 
seen in Fig. 9.7 (right). The solid red lines indicate applied cut. 

The missing mass distribution of two protons for events remaining after all previous 
cuts is still contaminated with background originating from the r] meson decay into pions 
and also direct two and three pion production (see Fig. 9.8). As the branching ratio is 
unknown for searched decay and also the total cross section for the direct three pion pro- 
duction remains unmeasured at the investigated beam energy it is not possible to estimate 
the absolute contribution from background. However, to remove as much background as 
possible we will restrict missing mass of two protons only to the narrow range around 
the rj meson mass from 544 MeV/c^ to 552 MeV/c^ (this range corresponds to the value 
of the FWHM of the resolution of missing mass distribution). The events which fulfill 
this condition will be considered further for the analysis. The simulated missing mass of 
two protons for background reactions and the signal is presented in Fig. |9.8| (left). The 
analogous distribution for the experimental data is presented in Fig. 9.8 (right). 

Finally, one can plot the invariant mass of particles in the final state: TT^e~^e~, and 
select the range corresponding the signal reaction. Figure |9.9| (left) presents simulated 
signal and background reactions remained after all previous selection conditions. One can 
see that the background reactions populate the lower invariant masses, and the signal is 
pronounced around the mass of the 7] meson. Therefore, to suppress the background we 
select invariant mass region from 555 MeV/c^ to 650 MeV/c^. This reduces the background 
nearly by 99%, decreasing the signal only by 53%. Later on this selection condition will 
be further optimized. The experimental distribution of corresponding invariant mass is 
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Figure 9.8: Distribution of the missing mass of two protons: (left) from simulations, and (right) 
measured data. The red sohd hne indicates selected region close to the j] mass in the range from 
544 MeV/c2 to 552 MeV/c^. 
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Figure 9.9: Invariant mass of particles identified as n'^e^e^ from: (left) simulations, and (right) 
measured data. The red solid line indicates the region of invariant masses from 555 MeV/c^ to 
650 MeV/c^ used to select candidates for the signal reaction. 



shown in Fig. 9.9 (right). The summary of all applied selection criteria with the number of 



experimental events remaining after the cuts and with the selection efficiency of the signal 
obtained from the simulation are shown in Tab. 19.21 

After application of all above described conditions to the experimental data number 
of event candidates for 7] — )• TT^e~^e~ decay is equal to N'^^^ = 10 it 3stat- However, the 
estimated number of N^^P may account not only for events corresponding to the searched 
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Cut 


Cut description 


£5 


Data 


PRESEL 


Preselection of two protons in FD, two 
oppositely charged particles in CD and 
two or more neutral clusters in SEC. 


13% 


7485153 


/M(77) 


Selection of two photons originating from 
the decay of vr*^ meson with invariant 
mass in the range from 90 - 180 MeV/c^ 


90% 


4952520 


iM{i+r) 


Restriction on invariant mass of two particles 
measured in MDC to be less than 40 MeV/c^ 


76% 


1154086 


SEC 


Particle identification in the 
Scintillating Electromagnetic Calorimeter. 


56% 


432813 


Split-off 


Suppression of events with the bremsstrahlung 
effect for charged particles. 


79% 


44668 


Conversion BP 


Suppression of external conversion of photons 
on the beam pipe. 


43% 


11617 


/M(77) 


Additional restriction on mass of vr" 
meson from 120 - 150 MeV/c^. 


69% 


8129 


MM(pp) 


Missing mass of two protons with the 
restriction of 544 - 552 MeV/c^. 


47% 


583 


IM(e+e-7rO) 


Invariant mass of final state particles with 
restriction from 555 - 650 MeV/c'^. 


47% 


10 


TOTAL 


Total efficiency for signal selection 
and final number of selected events. 


0.26% 


10 



Table 9.2: The list of conditions used for selecting the 77 e+e^7r" decay. Third column in- 
clude the efficiency reconstruction for the signal determined from the simulations. Fourth column 
indicates the number of experimental events left after application of subsequent cuts. 



decay r/ — )• vr^e^e" but it may also be due to the misidentification of the background 
reactions. In order to estimate the number of background events we have performed the 
same analysis of the simulated samples of background reactions (see Tab. 9.3). 

After application of selection criteria to the simulated background reactions, and af- 
ter taking into account the cross section and branching ratios values the only background 
channel left is pp — )• ppir^n^ — )• ppe^e~^^^ reaction which amounts to Nb = 13 it 4:stat- 
Therefore the number of reconstructed events in the measured data is consistent statisti- 
cally with the estimated amount of background. Thus we state that in the measured data 
sample of pp — t- pprj reaction we did not observe a signal from the searched 7] — t- 'ir^e~^e~ 
decay. In the next chapter we will estimate the upper limit for the branching ratio of 
investigated decay channel. 
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iNo. 


Reaction 


Number of 


Number of 


Number of 






generated 


reconstru- 


expected 






events 


cted events 


events 


1 


PP PPf] ppK^e^e~ — >■ ppe"'"e~77 (signal) 


5 X 10^ 


9208 


? 


2 


PP ~^ PPV ~^ ppTT'^TT^TT^ — >■ pp7r+7r~77 


20 X 10^ 








3 


PP PPV pp7r'^7T~Tr^ — >■ ^?p7r+7r~e"'"e~7 


1 X 10^ 








4 


PP PPV ~^ pp7r"^7r~7 


5 X 10^ 








5 


PP ~^ PPV ppe^e"^ 


20 X 10^ 


2 





6 


PP ~^ PPV ~^ pp7r°77 — >■ pp^^^^ 


3 X 10^ 








7 


PP ~^ PPV ~^ pp7r°77 — )■ ppe+e~777 


3 X 10^ 


13 





8 


jop — > ppry — )■ PP77 


6 X 10^ 








9 


— )• ppTT^TT~TT^ — )• pp7r+7r~77 


20 X 10^ 








10 


— )■ ppvT+vr^ 


58 X 10^ 








11 


PP — )■ ppTT^TT^ — >■ ppyj^J 


18 X 10^ 








12 


PP ~^ ppTT^TT^ — >■ ppe ''"e~777 


8 X 10^ 


7 


13 



Table 9.3: List of simulated reactions. Third column includes the number of initially generated 
events, and fourth column indicates the number of events left after all selection conditions. Last 
column indicates number of events expected to fulfill all selection criteria after taking into account 
a cross section and branching ratio. 
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10. Results for the ij Tr^e^e decay 

In this thesis we aim at estabhshing the unknown branching ratio for the rj — t- vr^e^e" 
decay which might violate charge conjugation invariance. However, as it was shown in 
previous section as a result of the conducted analysis of an experimental data for the pp — )• 
ppr] reaction we did not observe a statistically significant signal over the background. The 
final number of candidates events for the r] Tr^e'^e~ process selected from experimental 
events is equal to N^^p = 10±3stat, and the expected number of background events obtained 
from the simulations is equal to Nb = 13± Agtat- In such case, where the expected number 
of signal events is equal to zero, it is impossible to calculate the value of the branching 
ratio, and only an upper limit can be estimated. 

In order to calculate the upper limit for the branching ratio one has to know the 
number of expected: (i) signal, (ii) background, (iii) normalization channel events, (iv) 
selection efficiency, and (v) the branching ratio for the normalization channel. In case of 
this thesis for the normalization we choose the r] — )■ 7r'''7r~7r'^ decay. This channel has the 
same topology of the final state, as the signal reaction, thus the systematic effects should 
affect both processes nearly in the same way. Therefore, this will allow to cancel out many 
of the systematic uncertainties. 



10.1 The upper limit for the branching ratio 

BR{r] 7r°e+e-) 

The number of the signal events can be given by the following formula: 

Ns = L- ar,- BR^^^o^+e- ■ es, (10.1) 
and for the normalization channel: 

where L denotes the integrated luminosity, the cr^ is the total cross section for production of 
the r] meson, Ns and Nnorm stands for the number of events for the signal and normalization 
channel, eg and e„orm denotes the detection efficiency for the signal and normalization 
reaction, respectively. Therefore, the formula for the investigated branching ratio obtained 
from the two above given equations will take following form: 

D n ■ BR^^^+^-^0 ■ €norm , ^ 

B^V^wOe+e- = jij- ■ [W.6) 

99 



(10.2) 
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The branching ratio for the normahzation channel rj 



TT+vr TT^ is equal to BR{7] 



TT Vr TT 



(22.74 lb 0.34)% 12 . The signal for the searched decay rj — )• vr'^e+e was not 



observed in the investigated data sample, and therefore the equation (10.3) will take form 
of an inequality: 

BR^^^Oe+e- < jTr ' 

where the Ng^ denotes the number of expected signal events at a given confidence level. 

However, to test and optimize the last selection condition on the invariant mass of 
the 7r^e'^e~ we have varied this cut and for each case we counted the number of events 



for measured and simulated data. The result of these tests are presented in Tab. 10.1 
together with efficiency calculated for the signal. The cut discussed in Section 9.3 is shown 
as no. 2 in Tab. |10.1| One can see that by narrowing the condition on the invariant mass 



No. 


IM(7r0e+e-) 




£5 


jij-expected 


n— 


1. 


550 - 650 MeV/ 


'c^ 


12 


0.28% 


17 


1.64 


2. 


555 - 650 MeV/ 




10 


0.26% 


13 


1.96 


3. 


560 - 650 MeV/ 




7 


0.20% 


11 


1.82 


4. 


580 - 650 MeV/ 




3 


0.05% 


4 


1.25 



Table 10.1: Table showing the number of events obtained from the measured data and from 
simulations of the background channel pp — > ppn^Tr^ — > ppe+e^777 for different values of the cut 
on the invariant mass of the Tr^e~^e~ system. 



the number of events is decreasing both for experimental data and simulated background. 
Also the efficiency of the signal reconstruction decreases. As an criterion for choosing the 
best cut range we defined the significance parameter as: 



P 



N 



_e5 

expected ' 
B 



(10.5) 



The p parameter is increasing with the higher efficiency of the signal reconstruction and 
decreasing with number of background events. Thus as an optimum cut we have chosen 



cut number 2 listed in Tab. 10.1 for which the p coefficient was the highest. 

After the cut optimization one can estimate the upper limit for the searched branching 
ratio. The overall detection and reconstruction efficiencies for the signal and normalization 
channels were determined, based on the simulations and tliey are equal to: Cnorm — 

0.89% 

and €s = 0.26%, respectively. From the Tab. 10.1 one can see that we have measured 10 
events in experimental data and expect based on the simulations 13 background events. 
Thus, the value of the upper limit of the expected signal is equal to 



3.95 



123 



at the confidence level of 90%. The final result for the upper limit of the branching ratio 
equals to: 



< 



3.95 • 0.2274 • 0.89% 
82725 • 0.26% 



(10.6) 
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i?<4.0e+e- < 3.7 • 10-5 



,0,+,- ^ • xu (CL = 90%). 

Discussion of obtained result will follow in the next section. 



(10.7) 



10.2 Discussion 



Presented analysis for search of the 7] — )• ■K^e'^e" decay is the first from the WASA- 
at-COSY data where the rj meson was produced in proton-proton collisions. Obtained 
result is smaller than presently known upper limit given by the PDG group Il2| based on 



previous results 47 . This result, constitutes a next step in the search for rare decay of the 



rj meson by means of the WASA-at-COSY detector. Previously this decay was investigated 
by WASA-at-COSY in pd — t-^ Herj reaction where obtained value of the upper limit for 

10-5 



< 9 
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The result obtained in this 



the branching ratio was BR{r] 
work is compared to previous results and the theoretical prediction in Fig. |10.1 
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Figure 10.1: Comparison of obtained value of the upper limit for the branching ratio of the decay 
T] — )■ Tr'^e^e^ with previous measurements and the PDG fitted value. The red dashed line indicates 
the limit of the Standard Model predictions. 

In previous chapter, we have shown that analysis conducted in the framework of this 
thesis revealed no observe signal from the rj — )■ -K^e^e'^ decay. The experimentally obtained 
number of signal candidate events was the same as the number of estimated background 
events, with respect to the statistical uncertainty. Therefore, we determined an upper limit 
for the branching ratio for this decay. In case of small statistical signal there are many 



approaches which can be used for the estimation of the upper limits 1251. However, as a 



result of this thesis we will take the upper limit estimated in previous section according 
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to the prescription given by Feldman and Cousins 123 , assuming the Poisson distribution 
of the variables where the limit estimations on the physical quantities are based on the 
likelihood ratios. 



11. Summary and perspectives 



The main motivation of this dissertation was the test of the charge conjugation C in 
the strong and electromagnetic interactions. For this purpose we studied the rj — )■ tt^'7T~tt^ 
and r] ■K^e^e~ decays, where the r] meson was produced via pp — )• pprj reaction near the 
kinematical threshold. One of the goals was to determine the three asymmetries in the 
Dalitz plot population {Alr, Aq, As) for the strong isospin violating decay rj — >■ tt+tt^tt'^ 
in order to learn about the C invariance in different isospin states of the 7r+7r~7r° system. 
The second goal was to study the branching ratio for the r) — >■ Tr^e'^e" decay which can 
break the C symmetry in electromagnetic interactions, when occurring by forbidden first 
order transition. 

The measurement described in this thesis was performed in the Research Center Jiilich 
in Germany by means of the WASA-at-COSY detector installed at the Cooler Synchrotron 
COSY. The initial pp — t- pprj reaction was induced by the collision of a proton beam 
with the momentum of 2.14 GeV/c on a proton pellet target. Nucleons emerged from the 
interaction have been registered in the WASA forward detector, whereas the decay products 
of the ij meson were detected in the central detector. For the identification of the rj meson 
the missing and invariant mass techniques were applied. The background originating from 
the direct two pion production and other r) meson decays has been reduced to negligible 
level by applying the momentum and energy conservation lows, and by using the missing 
and invariant mass distributions, as well as by performing a kinematic fit method. For 
the T] —7- vr+TT^vr'^ decay the remaining physical background originating from the direct 
production of three pions via pp — >■ ppiT~^'K~iT^ was subtracted for each studied phase space 
interval separately. 

The three asymmetry parameters for the rj — t- Tr+vr^vr'^ decay, sensitive to the violation 
of the charge conjugation invariance were determined for the first time from the WASA-at- 
COSY data where the rj meson was produced in proton-proton reaction. The final sample 
contained about 10^ events which enabled to estimate the following values of the left-right, 
quadrant, and sextant asymmetry parameters: 

Alr = 0.0033 ± 0.0038,tat ± 0.0009,j/„ 

Aq = -0.0018 ± 0.0038,tot ± 0.0011,^/., 
As = 0.0006 ± 0.0039stat ± O.OOlOs^s- 

Established values of the asymmetry parameters arc consistent with zero within the range 
of the statistical and systematic uncertainty, which allows to conclude that the charge 
conjugation symmetry C is conserved in strong interactions on the level of the achieved 
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accuracy. Obtained result is also in agreement with previously measured values 26 26 39 



and the average of the Particle Data Group 39 



Furthermore the Dalitz plot distribution was analyzed. The Dalitz plot was represented 
in the X and Y variables proportional to the difference of the vr"'' and tt~ energies, and 
TT^ energy, respectively. The distribution was corrected for each bin for the geometrical 
acceptance of the WASA-at-COSY apparatus as well as, the reconstruction efficiency and 
was determined free from the background of the multi-pion production. The obtained 
density distribution was fitted with the phenomenological model. As a result of the fit 
we have obtained coefficients c = 0.05 it 0.10 which is standing in the odd powers of X 
variable, and it is sensitive for C violation. Obtained result is consistent with zero within 
estimated uncertainty. 

In second part of the thesis the investigation of the rj — )• 7r^e~^e~ decay is presented. 
Based on analyzed experimental data N^^P = 10 it Sgtat event candidates have been iden- 
tified from the two weeks of data taking. The background originating from the other rj 
meson decays was suppressed to a negligible level. The only remaining background which 
was left in the signal region is originating from the pp — )• ppTr^ir^ — ?• ppe^e~^^^ reaction. 
Next, we have estimated that the expected number of background events in the signal 
region should be equal to Nb = 13 it Agtat- Based on given number of measured events, 
number of obtained background, and the estimated uncertainty we conclude that no signal 
was observed and therefore we have established the upper limit for the branching ratio of 
the rj —7- TT^e'^e~ decay: 

Si2(r/^7r°e+e-) <3.7-10"^ (90% C.L.) 

Obtained result is slightly better than presently know value of this branching ratio from 



the fit of the PDG group 12 



Results presented in this thesis were based on the data sample of about 5-10''' r] mesons. 
As a perspective for further improvement of the results on both decays: r/ — t- vr+vr^vr*' and 
rj —7- 7r'^e"'"e~, it would be interesting to use full statistic sample, since the WASA-at-COSY 
currently has collected around 10^ rj mesons in proton-proton collisions, which is one of the 
world's largest data sample for the r] meson. Therefore, the investigations presented in this 
thesis will be continued. Such statistics should enable to lower the statistical uncertainties 
for the determination of the asymmetry parameters which were investigated in this thesis 
by a factor of ten. For the rj — )• vr^e'^e" decay it would be possible to lower the sensitivity 
of the branching ratio determination about an order of magnitude. 



"Anyone who has never made a mistake 
has never tried anything new." 

(Albert Einstein) 
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